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SCATTEEING IN DISCRETE RANDOM MEDIA 
WITH IMPLICATIONS TO PROPAGATION THROUGH RAIN 

Louis Joseph Ippolito, Jr. 


ABSTRACT 

An investigation of the multiple scattering effects on wave 
propagation through a volume of discrete scatterers was ac- 
complished. The mean field and intensity for a distribution 
of scatterers was developed using a discrete random media 
formulation, and second order series expansions for the 
mean field and total intensity derived for one -dimensional 
and three-dimensional configurations. Ihe volume distribu- 
tion results were shown to proceed directly from the one- 
dimensional results. 

The analyses demonstrated that either discrete or continuous 
techniques may be employed for the mean Reid and intensity 
expansions, as long as care is taken to insure non-overlapping 
scatterers in the formulation. 

The multiple scattering intensity expansion was compared to 
the classical ^single scattering' intensity and the classical 
result was found to represent only the first three terms in 
the total intensity expansion. The Foldy approximation to 
the mean field was applied to develop the coherent intensity, 
and was found to exactly represent all coherent terms of the 
total intensity. 

An incoherent intensity term, secular in L, in path length, 
was found which was not accounted for in the Foldy approxi- 
mation result or in the 'sin^e scattering' formulation. 

The study has demonstrated the feasibility of using discrete 
random media techniques for the determination of multiple 
scattering effects in propagation through a volume of discrete 
scatterers, and has provided some insight to the more gen- 
eral problem of multiple scattering in a rain volume. 
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SCATTEBING IN DISCRETE RANDOM MEDIA 
WITH IMPLICATIONS TO PROPAGATION THROUGH RAIN 

CHAPTER 1 
INTRODUCTION 


The influence of rain on the propagation of electromagnetic radiation has long 
been recognized as a potential problem for many applications in radio science 
and engineering. Investigators, beginning in the early 1940’s, noted that radar 
echoes at centimeter wavelengths were modified during rain conditions. Radio 
operators observed an increase in the level of noise and Interference when rain 
was present in the vicinity of the receiver. 

A comprehensive study by Ryde (1946) provided a theoretical basis for the meas- 
urement and prediction of meteorological effects on radio waves, and much sub- 
sequent work was based on his early efforts. 

The rain medium was recognized as a configuration of isolated dielectric spheres 
which absorb and scatter incident energy in accordance with classical field 
theory. The solution of the rain problem was based on studying the character- 
istics of the single rain drop and applying a simple summation of intensity to de- 
termine the bulk effecis of the rain volume on the propagating wave. 

Gunn & East (1954) ^plied the theory developed by Gustav Mie (1908) for the 
diffraction of an electromagnetic wave by a dielectric sphere and desciibed the 
scattering and absorption of rain by equivalent cross-sections for the spherical 
raindrops. Ihey found semi-empirical relationships between rainfall rate and 
attenuation which could be described by a simple power law dependence. 

Medhurst (1965) provided an extensive systematic analysis of rain attenuation 
for a wide range of frequencies and rainfall rates, and numerically calculated 
attenuation coefficients from the Ryde formulation. He determined the Mie co- 
efficients for spherical drops based on the measured drop size distributions of 
Laws & Parsons (1943). Medhurst compared his calculations to twelve sets of 
measured data and found that the measured attenuation tended to exceed the maxi- 
mum possible theoretical attenuation in most Instances. 

The advent of extensive microwave communications links in the late 1960’s pro- 
duced a large body of measured rain attenuation effects and studies extending 
from 2 to well above 100 GHz, Improved computer technology provided more 
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refined calculations for microwave attenuation calculations from the classical 
Mie theory, Setzer (1970), Chu & Hogg (1968), Satellite communicatioas Haks 
provided fuither information on rain attenuation at frequencies above 10 GHz. 

Two recent survey articles, by Crane (1971) and Hogg & Chu (1975), provide an 
extensive summary of the theory and measurements available on propagation 
factors influencing wave propagation through rain. 

The classical development of the theory of wave propagation through rain has 
progressed as basically a 'deterministic’ formulation, even though the phenomena 
involved are random in nature. The intensity of the wave propagating through 
the rain volume is assumed to decrease exponentially as it progresses through 
the precipitation. The quantities which describe the rain structure, such as drop 
size dist3?ibution, precipitation rate, refractive index, temperature, etc., are 
usually assumed as averaged quantities with deterministic characteristics. 

Since the decrease in intensity, or attenuation, of the wave propagating through 
the rain is deter mi ned as a simple summation of the power absorbed and scattered 
by each isolated rain drop, multiple scattering effects between drops are not ac" 
counted for directly. Because of this situation, the classical attenuation develop- 
ment has often been termed the 'single scattering’ approximation in the literature, 
for example Van de Hulst (1957), Medhurst (1965), Ishimaru & Lin (1972), even 
though no such explicit approximation appears in the theory. 

Some investigators have recognized the lack of multiple scattering considerations 
in the classical formulation, but disagree on its relative importance. Medhurst 
(p. 563) states that the neglect of multiple scattering effects along the path could 
be one source of error in the disagreement between his calculated values and 
experimental measurements. Hogg & Chu 1311) state that multiple scattering 
effects are impUcitly accounted for in the classical formulation since the incident 
wave at each layer in the rain volume talces into account attenuation by preceeding 
layers, therefore the bulk of the multiple scattering contributions are included. 
Crane (p, 179) states that the classical attenuation calculations strictly apply 
only to a nonscattering medium, and that neither adequate measurements nor an 
adequate theoretical treatment of multiple scattering effects have been made. 

Setzer (1970) applied the equation of radiative transfer given by Chandrasekhar 
(1960) and computed the first order Born scattering contribution for a line-of- 
site propagation link with a very narrow angular field of view. He concludes 
that, for a 150 mm/hr rainfall rate with a Laws & Parsons drop distribution, the 
multiple scattering term is completely negligible at 30 GHz, and of minor sig- 
nificance even at optical wavelengths. 

Crane (1970), using the same radiative transfer technique, but apparently with- 
out the assumption of a narrow field of view, presents curves (his Fig. 9) which 
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show that multiple scattering in a Laws & Parsons rain will be important at 
frequencies above 20 GHz for even moderate rain rates. The curves show, for 
example, that at the 150 mm/hr rain rate and for 30 GHz, multiple scattering ef- 
fects become significant at a path length of about 200 meters. 

The subject of multiple scattering has been treated extensively in the study of 
random media problems. In contrast with the ’deterministic' techniques em- 
ployed by classical rain attenuation methods, the study of random media re- 
quires a 'statistical' approach and 'averaged' results are of interest. The study 
of random media is approached by both continuous and discrete formulations. 
Frisch (1968) provides a comprehensive survey of wave propagation in random 
media for both the continuous and discrete scatterer cases. Major contributions 
in the continuous random media area include those of Keller (1962, 1963), 
Rosenbaum (1969), Kupiec, et al. (1969), among many others. 

A systematic procedure for the study of random media consisting of discrete 
scatterers was presented by Foldy (1945) and extended by Lax (1952), Bazer 
(1959), Lin & Ishimaaru (1971, 1974), and Twersky (1962-1967), among many 
others. 

The application of random media techniques to the rain scatter problem has not 
been very extensive, however. Lin & Ishimaru (1971) consider multiple scatter- 
ing effects on millimeter wave propagation in an isotropically scattering dis- 
crete medium. In a later paper Ishimaru & Lin (1972), the case of a narrow 
beam line-of-sight propagation path is assumed. Both studies present integral 
equations for the mean field and coherent intensity based on the Foldy approxi- 
mation formulation for forward scattering only, and solve equations after a 
number of simplifying assumptions which appear to reduce the problem to the 
classical, 'single scattering' formulation. The relative contributions of multiple 
scattering to the rain scatter attenuation are not discussed. 

The objective of this dissertation is an investigation of multiple scattering in 
propagaUon through a distribution of discrete point scatterers by a rigorous ap- 
plication of random media techniques. Chapter 2 reviews die continuous and 
discrete random media formulations important to the study. Chapter 3 establishes 
a treatment of propagation in one -dimensional media, with emphasis on the tech- 
niques applicable to the volume scatter problem. Chapter 4 extends the develop- 
ment to a volume distribution of discrete point scatterers, and develops expan- 
sion solutions for the mean field, coherent intensify, and incoherent intensify. 
Chapter 5 presents the general conclusions of the study of discrete scatterir^ 
in a random medium and discusses some implications of the results to the volume 
rain scattering problem. 
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CHAPTER 2 


WAVE PROPAGATION AND SCATTERING IN RANDOM MEDIA 


The subject of wave propagation in random media covers a vast range of topics 
and techniques, all sharing a common factor of statistical solutions to problems 
that axe non-deterministic in nature. As an electromagnetic wave propagates 
through a random medium, such as a turbulent atmosphere, fog, rain, ocean 
surface, or rough terrain, the randomness of the medium becomes the source 
of secondary waves, referred to as scattered waves or * scattering’. 

The scattering medium may be characterized in one of three ways; (a) a random 
continuum, where the randomness can be consider;-:d a continuous function of 
space and time, (b) discrete scattering, where the medium consists of a random 
distribution of discrete scatterers, and almowledge of the characteristics of a 
single scatterer is required, and (c) ron gh surface scattering, where a random 
surface function, rather than a volume fmction, is used. 

The study of rain scattering effects on ej.ectromagnetic propagation involves 
elements of both the random continuum and discrete scattering. This chapter 
reviews the basic formulation necessarj' for the developments which follow in 
Chapters 3 and 4. Most of the information in this ch^ter has been obtained 
from the published literature, and adapted, where necessary, for the particular 
notations and configurations of this study. The major references for the mate- 
rial summarized in this chapter are Frisch (1968), Foldy (1945), Twersky 
(1964), and Keller (1964). 

2-1. CONTINUOUS RANDOM MEDIA TECHMQUES 

A medium which is considered continuously random in position and time can be 
described statistically by a dielectric constant, , which is a random function 
of position and time, 


6j. = 6p(r, t) = n^(r, t) (2-1) 

where n^ is the index of refraction, also a random function. We malce the as- 
sumption that the time variations of the medium are slow compared to the time 
of wave propagation through the medium, which is a reasonable assumption for 
the rain scatter channel. Under this condition 

Cj. = 6j.(r) = n^Cr) (2-2) 

which implies that a propagating wave will maintain its frequency. 
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Consider Maxwell's equations for the medium, in complex vector notation 


V X E(r> = icou^HCr) (2-3a) 

V X HCr) = -iojeoe,(r)E(r) (2‘-3b) 

where a time dependence of is assumed. The electric field vector, for 

example, is available from 


e(r, t) = Re-|ECr) 

where Re-^ ^denotes 'real part of', and E(r) is a vector with complex 
components. 


(2-4) 


The vector wave equation for the medium is developed by combining Equations 
(2-3 a, b), 


V X V X E(r) - e^Cr) E(r) = 0 (2-5) 

with 


V • [e^(r)E(r)] = 0 


( 2 - 6 ) 


Equation (2-5) reduces to 


V2 E(r) + kl e^(r) E(r) - V 


V^) 


E(r) 


= 0 


(2-7) 


where 


The last term in Equation (2-7) can be neglected if the variations in over dis- 
tances comparable to the wavelength of the incident radiation are veiy small. In 
ne^ecting this term, polarization effects are also neglected since three uncoupled 
scalar wave equations result. 

The resulting vector wave equation is therefore 

E(r) + ^2-8) 

or in terms of the refractive index 

E(r) + n2(r)E(r) = 0 (2-9) 
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Each component of the vector wave equation above, ^^Cr) , satisfies Hie scalar 
wave equation 


^(r) + n2(r)^(r) = 0 


or 


[V^ + n2(r)3 ^^(r) = 0 (2-10) 

Consider a random medium imbedded is free space with a relative dielectric 
constant close to 1. The refractive index then has the form 

n^(r) = 1 + e u(r) (2-11) 

where the random function u(r) is defined by a mean <u(r)) and correlation 
(u(r)u(r' )> , and e is a small quantity. 

The wave equation then has the form 

+ k^[l + eu(r)]^^(r) = 0 

or 

[V^ + kjl ^(r) = -e k^ u(r) ^^(r) (2-12) 

This equation can be converted to an integral equation by use of the Green’s 
function method. The solution is assumed as the sum of an incident wave which 
obeys Equation (2-10) and a scattered wave. The associated boundary conditions 
require that the scattered wave be composed of outgoing waves at infinity, and 
that ^(r) is continuous and vanishes at infinity. With these boundary conditions, 
the wave equation converts to 

>k(r) = %(r) - 6 kj ; G(r, r') u(r') 'ff(r') dr' (2-13) 

where G(r, r') is the free space Green’s function 

_gik |r - r'l 

G(r,r')= (2-14) 

43r!r-r'l 


which is a solution to the wave equation 

(V2 + k^)G(r,r') = 5(r - r') 


(2-15) 
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The inoident wave ^o(r) is a solutLon to the homogeneous wave equation 

(V^ + k2)^^(r) = 0 (2-16) 

Most of the field analyses in this study will begin with the integral form of the 
wave equation discussed above. 

Born Approximation 


An approximate solution for the integral form of the wave equation for scattering 
in random media is available by assuming the solution is of the form 


oo 


n = 0 


(2-17) 


The first iteration solution, for n = 1, from Equation (2-13), is 

^^(r) = ^o(r) - ek^ / G(r, r') u(r') ^o(r') dr' + 0(^2) (2-18) 

where O(e-) indicates additional terms of order . 

The second iteration assumes Equation (2-18) for a solution to Equation (2-15), 

i.e. 


’^^(r) = ^^(r) - ekl f G(r, r') u(r')[^o(r') 

- 6k^ / G(r',r")u(r")^o(r")dr"l dr' + 0(6^) 


(2-19) 


or 


^(r) = %(r) - ek^ / G(r, r') u(r') ^^(r') dr' 

+ k^ / G(r, r') G(r', r") u(r') u(r") ^^ir”) dr' dr" + 0(e^ ) 


( 2 - 20 ) 


Additional terms of the Born series can be ^plied to obtain higher order ex- 
pansions. For this study expansions to O(e^) will be utilized. 


The expectation, or average, value of ^(r) is available from Equation (2-20) 


<*(r)) = 'F„(r) - ek^ / G(r, r') <u(r')> 'l^o(r') dr' 

+ f G(r, r') G(r',r") (u(r') u(r"))^o(r") dr' dr" + OCe^) 


(2-20a) 
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where <u(r')> and <u(r' ) u(r")> are the first and second order statistiCB which 
define the function u(r). 


2-2. DISCRETE RANDOM MEDIA TECHNIQUES 

Consider a distribution of identical isotropic point scatterers in the presence of 
an incident wave The scatterers are located at rj . The scattering prop- 

erties of the individual scatterer are characterized by a scattering function, gj , 
which is a deterministic function determined from the physical characteristics 
of the individual scatterer. The distribution of the scatterers is given by the 
distribution function f(r). We assume that g-, f(r), and are known. We 

require the average of the wave function, <w(r)) and other related quantities. 

For a particular configuration of scatterers, the wave function at the point r 
can be represented as (Foldy, 1945, p. 110) 

❖(r) = ^ gj C(r, rj) (2-21) 

j 

This equation represents the field as the sum of the incident wave and spherical 
waves diverging from each of the scatterers. The external field acting on the 
j th scatterer is then 


^(ij) = ^ gj-\E^(r^,)G(r,r.,) (2-22) 

j' 

j' j 


Equations (2-21) and (2-22) represent the fundamental equations of multiple 
scattering. 


The first order approximation to the solution of the above equations assumes 
that the external field on the j' th scatterer, ’l^(rj»)» is tJis incident field at that 
scatterer, i. e. 


'ir(r.,) ^ (2-23) 

The field at r is again represented as the sum of the incident wave and spherical 
waves diverging from each of the isotropic scatterers, i.e. 


^(r) = 4-oW + Z) G(r,rp (2-2S) 

3 
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where iSfCrj) is the external field acting on the j th scatterer, and G(r, rj ) is the 
Green's function representing the spherical wave of the j th scatterer. 

But now, in the first order approximafion, ^(rj) is the sum of the incident wave 
at Tj plus the sum of the fields produced by the other scatterers excited by the 
incident wave, i. e. 


^(rp = ^r^Crp + gj'G(r-,r-,)^oC^j') ( 2 - 26 ) 

j' 

j'^j 

where G(ij , rj/) is the field produced at rj from the scatterer at tj, excited by 
the incident wave ❖^(rj,). 

Substituting Equation (2-26) into Equation (2-25) 


= %(r) + 2 ^G(r,rj) 

j 


Sj.Gfrj.rjO'I'ofrjO 

i' 


or 


t 

2 gj G(r, i-j) + J] ^ gj. G(r, r^) Gfrj, r^,) 

j j j' (2-27) 

where the prime in the summation indicates j Equation (2-27) represents 

the fundamental equation for the evaluation of multiple scattering by discrete 
isotropic scatterers, and will be the starting point for the discrete formulation 
expansion solutions in this study. 

Foldy Approximation 


The direct method of solution for the discrete multiple scattering field, repre- 
sented by Equation (2-27), consists of solving the set of simultaneous linear 
equations for the ^(rj), given by Equation (2-26), and substituting these into 
Equation (2-27) to find ^(r) as a function of position and the parameters of the 
scatterers. This is not a feasible method. The problem is even more severe 
when <^f(r)) is desired. In that case, integrations must be carried out over all 
the simultaneous equations to obtain ('^(r)) . 

This problem was studied by Eoldy (1945), and resulted in the development of a 
procedure for determining the mean field integral equation by applying an ap- 
propriate approximation for the external field acting on each scatterer. This 
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important development has been used extensively to study multiple scattering 
effects, and serves as the basis for most present treatments of discrete scatter- 
ing problems, Twershy (1964), Ishimaru & Ldn (1972), Lax (1951), Baser (1959), 
among others. 

The detailed development of the Foldy approximation will be applied in this study 
for the one-dimensional case, Section 3-3, and for the volume scattering case, 
Section 4-2. The Foldy approximation ejqjansion for <^(r)) and K^(r)>l^ mil be 
compared mth directly evaluated expansions for the slab and volume scatter 
cases, and the limits of applicability to the rain scatter channel evaluated. 
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CHAPTER 3 


PROPAGATION IN ONE-DIMENSIONAL MEDIA 


As a prelude to the direct evaluation of multiple scattering in a volume of scat- 
terers, we begin by first establishing a rigorous treatment of the problem of the 
multiple scattering of waves by distributions of point scatterers. In ibe next 
chapter, similar techniques will be applied and extended to the three-dimensional 
volume scatter problem. 

A one-dimensional slab of Poisson distributed point scatterers will be considered, 
and expansions of the mean Held and intensify will be developed, utilizing con- 
tinuous and discrete techniques. The particular areas of interest to be investi- 
gated include, (a) the relative advantages of discrete versus continuous tech- 
niques for point scatterer problems, (b) the applicability and range of validity 
of the Poldy approximation in the mean and intensity expansions, (c) the deter- 
mination of a 'weighting factor' for the point scatterer which represents a di- 
electric finite scatterer in the limit as its width approaches zero, and, finally, 

(d) the implications of the one-dimensional results and their applicability to the 
solution of the rain scatter problem. 

Many of the techniques and methods of evaluation established in this chapter will 
be found to be directly applicable, with some extensions, to the three-dimensional 
volume scattering problem, hence the results here will provide a useful guide in 
the method of analysis for this study. 


3-1. POISSON IMPULSE REPRESENTATION FOR DISCRETE SCATTERERS 


The point scatterer will be utilized in this study to represent a physical scatterer, 
thus it is important, at least in the mathematical sense, for its properties to 
reasonably represent the characteristic properties of the scatterer. It will be 
useful in the statistical evaluation of the moments of a configuration of point 
scatterers to consider the scatterers as a uniform distribution of Poisson im- 
pulses, with an appropriate 'weighting factor' which represents the physical 
drop characteristics, such as radius and dielectric constant. In this way the 
point scatterer formalism is maintained for ease of mathematical evaluation, 
yet the physical properties of the scatterer are maintained and available for 
practical interpretation of the results obtained. We begin by considering a 
single point scatterer and developing the 'weighting factor' which represents a 
dielectric scatterer within specified bounds of size and dielectric constant. 
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(a) Single Delta- Function Scatterer 


Consider a single impulse scatterer located at in free space, as in Figure 
3-1 (a), on the x axis, with an incident wave propagating from the negative x 
direction of 


^j(x) = (3-1) 

where an e-^^^^time dependence is assumed. The field "Fix) obeys the reduced 
one-dimensional wave equation 


— r ^(x) + n^(x) ^FCx) = 0 

dx^ 


(3-2) 



X 


(a) 



(bj 

Figure 3-1. One-dimensional Scatterer Configurations 
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where n(x) is the index of refraction and k is the free space wave mimber. The 


scatterer at results in 


n^(x) =1+0 6(xi) 


where 5(X|) is the Dirac delta function and a is a complex 'weighting factor', to 
be determined later. The wave equation is therefore 


'P(x) + [k^ + k^o 6(Xj)3 = 0 


— r + k^ ^(x) = -k^o 6CxO 


The above equation converts into the integral equation 

^f(x) = ^j(x) - k^o J Gjj(x, x') 6(x' - Xi)>k(x')dx' 
where G (x, x') is the free space Green's function, i.e. 


Jk|x -x'l 


Gq(x,x') = 


The delta function evaluates the integral at x' = Xj 

giklx - xil 

^(x) = ^j(x) - k^a ^(xD 

2ik 


At X = X; 


ikx* 

^(x.) = e ''''' - — ^{xD 
‘ 2i ^ 


, ‘ ^ ikXi 


(3-S) 


Substitute ^(k.) back into Equation (3-7) 


, ik|x-Xj! 
ko e ’ 


^(x) = ^j(x) - — 

‘ 2i 


ikX; 


ka 

1 + — 
2i 


(3-9) 


For X > X. 


ka 


^(x) = ^^j(x) - — 


,ikx 


2i 


ka 
1 + — 
2i 


or 


^(x) = 


I 


ka 

1 -i — 
2 


likx 


(3-10) 


The term in brackets is seen to be the transmission coefficient Tg, of the delta 
function seatterer, i.e. 


Ts = 


1 


1 -i 


ka 


(3-11) 


(b) Evaluation of Weighting Factor from Transmission Coefficient 


In this section the transmission coefficient for a slab of random media will be 
compared with the delta function transmission coefficient and a weighting factor, 
a, developed. 


The transmission coefficient, T, for a dielectric slab of width L was developed 
in Appendix A, Equation (A-12) 


T = 


4kKe‘^L 


(k + K)2 - (k - K)2 e^2KL 


(3-12) 


where K is the propagation constant in the slab and k is the free space value (see 
Fig. 3-1 (b) for the slab configuration). Let be the complex relative dielectric 
constant of the slab, then 


K = k\/ir 


(3-13) 
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1 


r 


1 


We impose the following conditions on e, and L: 

(1) 62 is > 1, but 

(2) L is small, such tliat 

(3) kV^ L KL < 1. 

E3q)anding the denominator terms in T 


Ck + K)2 = k2(l + ^ e 




Ck - Kf 


= k^a - Vi;)2 = k2 e 




From condition (1) above 

{• *^)‘ ■ ■ 

and 


(‘ ■ v^) ~ ‘ ' vs; 


From condition (3) 

giKL ^ j + L 

and 

gi2KL ~ 1 + i2kv^ L 

Therefore, under the above conditions, T will be 


4k^ (1 +ikVi^ L) 

k^ ^2 

1 

V vs;y 

l-l 


|(1 +i2kVi^ L) 


(3-14) 
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4(1 + ikVeT L) 

T = ^ 

4 + (l -V^)i2kVi^ L 


(3-15) 


With conditions (1) and (3) employed above, T reduces to 


4 

T = 

4 - i2k 62 L 


and finally 


1 

T ^ 

k e,, L 
1 -i — ^ 
2 


Comparing this result with Tg, Equation (3-11) 

1 1 
o 

ka ke- L 

1-i — i-i — ^ 
2 2 


we see that 


a = £3 L (3-18) 


Thus, with the weighting factor £2L, a delta function scatterer can represent a 
physical dielectric scatterer of dielectric constant gj and width L, under the 
limiting condition 


kVe^ L 1 (3-19) 


With k replaced by 2w/'K, this condition requires that 


2irL 1 

< ^ 1 (3-20) 

X v^2 


(3-16) 


(3-17) 


which is seen to be the Rayleigh scattering criteria, i. e. the size of the scatterer 
must be much smaller than the incident wavelength. 


(c ) Mean Field. - Continuous Formulation 


In this section we consider the propagation of a plane wave incident on a slab of 
random medium, where the medium is described by a continuous function of 
position and time. This description, where the medium can be considered a 
"random continuum, " differs from the discrete scatterer description, which will 
be developed in the next section. 

Figure 3-2 shows the configuration under consideration; a slab of width L, with 
complex refractive index n(x), and a plane wave incident from the negative x 
direction. The refractive index of the slab can, in general, be a random funcli:>n 
of position and time, n(x, t), but we assume that the time fluctuations of the me- 
dium are slow compared to the time of wave propagation, therefore the refractive 
index can be considered a function of position only. 

The scalar field ^'(x) satisfies the reduced one- dimensional wave equation 

d2 

— ^(x) + k2 n2(x)’^(x) = 0 (3-21) 

dx^ 

We assume n^ (x) in the slab is 

n^(x) = 1 + eu(x) 0 < X < L (3-22) 



Figure 3-2. Plane Wave Incident on a One-dimensional Slab 

of Random Medium 
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where u(x) is a random function of x, defined by a mean<u(x)>, and correlation 
<u(x) u(xf)). Let the function u(x) be a Poisson impulse process of the form 

N 

u(x) = ^ a 6(x - Xj) (3-22a) 

j= 1 

Then the mean and correlation are (P^oulis, 1965) 


<u(x)) = crX 

<u(x) u(x')> = + a-X fi(x - x') 

where X = N/L, the uniform number density. 

The wave equation is then 





1?(X) = -k^ 6 u(x) ^(x) 


(3-23) 

(3-24) 


(3-25) 


The above wave equation and its associated boundary conditions convert into the 
integral equation 


^(x) = >^kQ(x) - 6 


i 


L ^ik'ix “ x'i 


2ik 


u(x') ^'(x') dx' 


(3-26) 


where ^ (x) is the incident wave, a solution to the homogeneous wave equation 


dx^ 


+ k^ 


^(x) = 0 


(3-27) 


To obtain the Born expansion solution, let ^(x) be of the form 




(3-28) 


n - 0 


The first term of the expansion gives 

.L „ik|x - x'l 


’Sr(x) = - k^ e 


X 


2ik 


u(x') dx' -t O(e^) (3-29) 
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Substitute ^(x) from Equation (3-29) into tbe integral of Equation (3-26) 


e 


I 


L Jk|x-x'| 


2ik 


u(x') 


s.ikx' 


- 


L 


L piklx'-x' 


2ik 


u(x") dx" 


dx' + O(e^) 


or 


^(x) = - k^ e I 

J n 


L pik|x - x'l 


2ik 


u(x')e^^'dx' 


+ k'^ 


i 


L pik|x-x'l 


2ik 


u(x') f 

J n 


L giklx -x”| 


2ik 


(3-30) 

u(x")e^*^" dx^'dx' + OCe^) 


Combining the exponent terms 

•L eik(lx - x') 


^f(x) = -k^ e 


/ 

•'o 


2ik 


u(x^) dx' 


(3-31) 


+ k'^ 


'■// 

*/ o «'o 


L /-L pik(|x - x'l + lx' - x"l + x") 


~4k^ 


u(x') u(x") dx' dx" + 0(e3) 


Equation (3-31) gives ’Jf(x) to O(e^). 


Additional terms of the Born series can be applied to obtain higher order ex- 
pansions, i. e Equation (3-31) is substituted back into the integral of Equation 
(3-26), etc. , to obtain ^r(x) to O(e^). For this study ^(x) to O(e^) will be 
utilized. Taking the average of Equation (3-31), we have 


<^r(x)> = - k^ e 


i 


L pik(lx - x'l + x') 


2ik 


<u(x')> dx' 


(3-32) 


+ k'*^ 


L pik(Ix - x'l + lx' “ x"l + x") 


// 

“'o “'O 


(u(x') u(x")> dx' dx" 


-4k^ 
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Substitute in the mean and correlation of u(x) from Equations (3-23) and (3-24) 





gik(lx -x'l + x') 


4 



gik(|x - x'l + lx' 


X I + X ) jjjjf jjjj" 


(3-33) 


+ X 



gik(|x - x'; 


|x' - x"| + x") gf'x' - x") dx' dx" 


The delta function in the third integral evaluates the integral at x' — x", i.e. the 
third integral term is 


4 



g,ik(lx - x’l + x') 


The resulting mean field is therefore 


ke 

<^^(x)> = aX 

2i 



eik(|x-x'| + x') 


4 



eik(lx-x'l + x') 


(3-34) 


(3-35) 


+ 



gik(|x - x'l + |x' - x"| + x") 


The above equation expresses the mean field inside and outside the slab, to 
O(e^), and was derived using the classical Born expansion techmque for an as- 
sumed continuous random medium* 


(d) Mean Eield - Discrete Formulation 

We consider here the same one- dimensional slab problem treated in Section (c), 
under the assumption that the random medium consists of a uniform distribution 


22 









of N discrete scatterers whose position Xj is defined by the Poisson impulse 
process 


N 


a “ Xj) 0 < X < L 

1 


(3-36) 


where a is a complex weighting function defined by the characteristics of the 
scatterer. 

We proceed to find the solution for the field ^(x) for any value of x. The field 
satisfies the reduced homogeneous wave equation outside the slab 


2 

— - + 
.dx^ 


^n(x) = 0 


X < 0 
X > L 


(3-37) 


and 

>ifQ(x) = 

In the slab, the wave equation for the field is 


(3-38) 


— - ^(x) + 
dx'^ 


N 


i + ^ 5(X - Xj) 

j = 1 


^(x) = 0 


(3-39) 


or 


'-£l 

_dx2 


+ 


N 


^f(x) = “k^o ^ 6(x - Xj) i^(x) 
j= 1 


(3-40) 


The above wave equation and its associated boundary conditions convert to the 
integral equation 


N rh gik|x - x'l 

^(x) = ^fp,(x) - k^ff I ^ 6(x' -Xj)^(x')dx' 

j = 1 0 


(3-41) 


The delta function evaluates the integral at x' = x. 


kff N 

^(x) = ^^q(x) + i — ^ gik|x-Xj| ^(x.) 

2 j=i 


(3-42) 
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This equation represents the field at x as the sui^of the incident field in the ab- 
sence of the scatterers, and waves of the form e* with an external excit- 

ing field '^{Xj ) erainating from each scatterer. 

The external field at the j th scatterer is the sum of the incident wave at x. plus 
the sum of the fields produced by the other N-1 scatterers, i.e. 


ikx; 


ka 


N 


Jk|Xj - xjl ikxj 


^f(Xj) = + i — ^ e ^ ‘ e 

^ 1=1 
j 


(3-43) 


ilc jx * “ X 1 1 

where e ^ ' is the field produced at xj from the scatterer at x^ excited by 

the incident wave e^^^^ 

Substituting '4^(Xj) from Equation (3-43) into Equation (3-42), 




2 j=I 


ikXs , , iklx;-x,| -,y \ 

® ^ I — 2^ ^ 

2 I 


N 


or 


IffT ^ *1 I \ rr^ ^ ^ ^ i r i ^ 

^(x) = ^ ^ ^ “ ""j' ’""j ■ 

^ j = 1 ^ j = 1 I = 1 (3_44) 

where the ' denotes 1 7 ^ j . 

It is of interest at this point to compare this result with a Born e^ansion solu- 
tion of 4^(x), Equation (3-42). Assume 


*(x) = X) o" 


( 3 - 45 ) 


n = 0 


The first term of the expansion, from Equation (3-42), gives 


4-(x) = e>kx + i 2 -’=i' e“=*i -i- OCa^) 


2 i = I 


(3-46) 
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Substituting from Equation (3-46) into the summation term in Equation 
(3-42) 


4^(x) = + 



ik|x - Xii 
B J 


ikxf 
e J 


+ 




0(o3) 

(3-47) 


or 


^(x) = + i — 

2 j=i 


^ik(lx - xjl + Xjj 


k2 (j2 

4 


N N 
j = 1 1=1 


^ik(lX-Xj|+ IXj-Xil + Xj) ^ 


(3-48) 


Comparing Equation (3-48) with Equation (3-44), flie results are seen to be iden- 
tical, except for the prime in the summation of Equation (3-44), i.e. 1 j. Thus 
the Born expansion represents the discrete scatterer field if 1 7 ^ j, i. e. if the 
scatterers are not allowed to overlap. The significance of this point will be dis- 
cussed in more detail later in this section. 


Proceeding to the average of Equation (3-44), we have 


C^f(x)) = + 



ik(lx -Xjl 


4 



^ik(|x - Xjl + ixj - Xjl + Xi)v 


(3-49) 


The evaluation of the mean terms above is accomplished in Appendix B, Section 
(a). From Equation (B-9) 



^ik(|x - Xjl + 




^ik(|x - x'l + x') 


(3-50) 
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and from Equation (B-18) 


(tt ^ik(|x - Xjl + |Xj - xil + xA „ 

\i=ll=i / 


rr-* 

-'o *^0 


(3-60) 

ik(!x - x'l + lx' - x"| + x") 


where 


X 


N 

L 


The resulting mean field is then 



gik(|x - x'l + x') jjjj' 


gik(!x - x'l + !x' - x"| + x") dx" 


(3-61) 


Equation (3-61) represents the mean field for a set of uniformly distributed non- 
overlapping point scatterers whose position is described by the Poisson impulse 
process of Equation (3-36). 

Comparing the discrete formulation result above with the. mean derived from the 
continuous random medium formulation. Equation (3-3S) with e = 1, the results 
are seen to be identical except for one term which appears in the continuous 
result, 


k2 

4 



gik(|x - x'l + x') jjx' 


This term is produced for the condition where the N scatterers are allowed to 
overlap, i.e. when x' = x", which is allowed in the continuous formulation. The 
term is produced from the delta function term in the correlation for u(s), Equa- 
tion (3-24), and results in the third integral of Equation (3-33), which is evaluated 
by the delta function at x' = x", Equation (3-34). 
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Note also that the Born expansion for the discrete formulation. Equation (3-48), 
allows overlapping seatterers and would also produce the term above. 

The comparison of the continuous and discrete formulations has brought out a 
very important factor relating to the condition of overlapping seatterers. For 
application to the physical rain scatter channel, the overlapping scatterer is not 
acceptable, and the inclusion of the overlapping scatterer term in the mean field 
expansion would lead to incorrect results. 

Therefore, if the continuous formulation is used, the correlation function for 
u(x) must be modified to 


<u(x) u(x')> = (3_62) 

which can be shown, (Papoulis, 1965, p. 285), to represent the Poisson process 
with non-overlapping intervals. The continuous formulation mean field would 
then be identical to Equation (3-61), the discrete result. 

One final point, regarding the discrete formulation utilizing the Born expansion; 
since overlapping seatterers are not allowed, the Born result Equation (3-48) 
must be modified with a prime in the double summation to produce the desired 
field. 

In summary, as long as care is taken in the formulation of the solution to Insure 
non-overlapping seatterers, either discrete or continuous techniques can be used 
for the development of the mean field and its associated parameters. 


3-2. SERIES EXPANSION SOLUTIONS FOR THE ONE -DIMENSIONAL SLAB 

The mean field was developed in Section 3-l(d), Equation (3-61), for a slab of 
width L consisting of N uniformly distributed non-overlapping point seatterers 
whose position is described by the Poisson impulse process of Equation (3-36), 
In this section we proceed to evaluate expansions to 0(ff^) for the mean field 
<^(x)>, coherent intensity Of (x)> <’f*(x)>, and total intensity <^(x) ‘f*(x)>, for 
the region x > L, 

(a) Mean Field 

Consider the Integral form of the mean field. Equation (3-61), We desire a 
solution for x > L, therefore with x > xj, xj 

(lx - Xjl + Xj) “ x ~ Xj + Xj = X 

(lx - Xjl + Ixj - Xjl + X| = X + IXj - Xjl + Xj - Xj 


27 


Then 


lc(j 

(■4r(x)> = + i — X 

2> 




(3-63) 




n L 

I 


,ik(lx' -x"l + x" - X ) dx' dx" 


Note that 


r 


dx' = L 


The double Integral is evaluated as Ij in Appendix C, Equation (C-5) 

:,ik2 L L 


Ice 

+ i — X^ 

2 4 


i 

+ 


\_2 4k^ 4k^ 2ikJ 


or 


<^(x)> = 


kff k"a2 Jl? 1 L 

l+i^XL-— + 4k, 


(3-64) 


(b) Coherent Intensity 

The complex conjugate of the mean, from Equation (3-64) above, is 


<4r*(x)> = 


kcr=>= k2(a^)* /L^ 1 L \ 

1 . i XL — X (-^ +—2 2ik) 


(3-65) 


The coherent intensity is found as 


ka 


<^(x)> <^*(x)) = K^(x)>p = 1 + i — X L - 

kcr* 


k2 JV- 1 


,ik 2 L 


X^ — + 


- i XL - 


2 4 \ 2 4k^ 4k^ 

k^(o^)* 1 


L> 

2ik> 


X" 


2 

k^ (j o* 


1 2 4k^ 


-ik 2 L L \ 


4k^ 2ik/ 


+ 


X2 + 0(a3) 
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or 


k (7 a* _ - 

= I + i — X L(cr - a*) + 

2 4 

k^a2 ^ /l 2 ^ 1 ei^2L l\ 

4 \2 4k2 4fc2 2ik/ 

fc2(t;^)* -/L2 1 g-ik2L L \ 

X^f — + + ) 

4 \2 4fc2 4k2 2ik/ 


(3-66) 


(c) Total Intensity 

The total intensity, (S&'Cx) ’^*(x)> , Is developed from the discrete formulation 
for ^(x), Equation (3-44). Taking the complex conjugate 


’I''*(x) = 



-ik([x - Xj'l + XjO 


k2(ff2)* 

4 



-ik(lx - Xj'l + |Xj< - Xj'l + Xj') 

0 


(3-67) 


Then 

^00 ^*(x) 


1 + i — e' 
2 


,-ikx ^ Jk(|x-xj| + xj) 


k2 (j2 


j = 1 
N N 


E E 


' ik(|x -Xjl + IXj -Xjl + Xj) 


j = 1 1=1 


y g-ik(|x-Xj'| + Xj') 

^ik(|x -Xji + Xj) y ^-ik(|x-Xj'| + Xj') 

^ ^ r=i 


k^ a a* 

N 

4 

j = 1 

k2 j* 

Jkx 


(3-68) 


j' = i r = i 
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The expectation value of Equation (3-68) is 


Wx) ^ ‘(x)) = 1 + i " ’‘l' 


4 

k2(a2)=P 


a”ikx 


' N N \ 

S - Xjl + Ixj - Xjl + Xj)\ 

J=I1=1 / 


(3-69) 


N N 


j = 1 1 = 1 


^ k^ (T ^ ^ gik(lx-xj| + xj-|x-xj-l-xj0^ 

J= I j'= 1 

The mean potions above are evaluated in Appendix B, Equation (B-9) for the 
first two terms. Equation (B-IS) for the next two terms, and Equation (B-17) 
for the last term. Applying these results 

Ico* /* ^ 

Wx)’|r=i=(x)) = (|iir(x)|)2 = 1 +i_Xe-ikx I gikCix - + x') dx' 


kg'*' 

- i Xc' 

2 


“X' 


Q-ik([x -x'l + x') 


k2 g2 
4 

k^(g^)'*' 

4 

a a* 


>? e-^kx r r"* 
0 Jq 


gik(|x - X i + |x - X '1 + x") Jjj_r 


J^2gikx I i e-ik(|x-x'!+ |x'-x"| + x")dx'dx" 


Li 

f-rr ' 


(3-70) 


ik([x - X I + x' - |x - x^'l - x") 


+ X 


j -x'l^x'-lx -x'|-x')^x^ 
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1: 


4- 






where X = N/L, Solving the integrals in Equation (3-70), and recalling that 

X > x', X > x" 


we get 


ka ka* 

<l^f(x)P> = 1 + i — X L - i X L 

2 2 


a f f dx' dx" 


If 

0 *^o 


(3-71) 


k^Ca^)* - i~ i~ .."I,.." 


^ f C 

Jo Jo 


^-ik(|x' - x"l + x” - X ) dx" 


k^ tJ £7* 


+ 


X^ 1? + 


k^ 


a a^' 


XL 


4 4 

The first double integral is evaluated in Appendix C, Equation (C-5) 


L = — + 


1 L\ 

4k^ 4k^ 2ik/ 


(3-72) 


The second integral is the complex conjugate of above, evaluated in Appendix 
C, Equation (C-7) 


/L2 1 Q-i2kL L' 

I = . — + + — 

^ V 2 4k^ 4k^ 2ik, 


(3-73) 


Therefore 


k jv- ! ei2'='- L\ 

2x = 1 4. i— \ TY/t _ nr*’» X^ — + T 

4 \2 4fc2 4k2 2ik/ 


<I^(x)|2> = 1 + i — XL(o - 0 =**) 
2 


(3-74) 

^ _L . . LV — i-^ - — L 

4 \ 2 4k^ 4k^ 2ik/ 4 4 
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or 


- k or (7* 

<|^(x)|2> = 1 + i - A L(cr - o*) + X L 

2 4 




k^ 

+ — 
4 


a 


+ (p^ )*\ + (tj^ )'■' 


L2 


4k2 


£1 ei^kL <ZZ e-i2kL 


4fc2 


4k^ 


2ik 


(3-75) 


Witli 


c; == Oj. + I ffj 

ff “ cr^' = i2ffj 

= 2(or^ - ff?) 
(7 CF* = or^ + (7? 

r 1 

(7^ - (CT“)* = i 4 cTj, CTj 


Thus Equation (3-74) is further reduced to 


k^ k^ 

(|■^r(x)P> - 1 “kXLCTj+ — X(ct^ + ff?) L + — X^ 

14^1 4 


io^-<yf) o a- 

+ r~ cos 2 k L -t- — — sin 2 k L 


(a^-ah 2a a. 

2k^ k 

(3-76) 


2fff V 


2k^ 


4k^ 


The above result is the total Intensity in the region x > L for the slab of width 
L. For L = 0, <l'4r(x)P> = 1, as expected. We further note that for a = ia^ 
(pure complex) 


k^ k^ 

<|^(x)j2) = 1 - kXLa. + — Xa? L + — X^ 

1 4 1 4 


And for cr = (real) 


2a? (1 “ cos 2 k L) 

^ 2k^ 




(3-77) 


k^ k^ 

<|^r(x)|2> =1+— XffH + — X2 
4 4 


L (1 - cos 2 k L) 

„ 2k2 


(3-78) 
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If the total intensity expansion, Equation (3-Y4), is compared with the coherent 
intensity expansion, Equation (3-66), all terms are the same except for 

Q a* 

-i- XL 

4 

which is present in the total intensity expression. This term represents the 
first term in the incoherent intensity expansion, i, e. 

(I’S^(X)P) - K^f(x))P = — — XL -{- 0(ff3) (3-79) 

4 

and is seen to be of O(ff^L), 


3-3. APPLICATION OP THE FOLDY APPROXIMATION TO THE ONE- 
DIMENSIONAL SLAB 

The problem of multiple scattering of particles by a random distribution, treated 
by Foldy ( 1945), resulted in the development of a procedure for determining the 
mean field integral equation by applying an appropriate approximation for the ex- 
ternal field acting on each scatterer. This important technique, later utilized 
by Lax (1951), Twersky (1964, 1967), Ishimaru & Lin (1972), and others, is 
reviewed in Section 2-2, 

The Foldy approximation is applied here to the one-dimensional slab problem, 
the mean field <^(x)> developed, and the resulting expansions compared to the 
expansions previously obtained. 

Consider the one-dimensional wave equation for the slab of width L, consisting 
of a uniform distribution of N scatterers, given by Equation (3-40) 

N 

\f(x) = -k^ a 5(x - Xj)^f(x) (3-80) 

j= 1 

The above wave equation plus the associated boundary conditions produce the 
integral equation 



^f(x) = 



giklx - x'l 

5(x' - xO 4^(x') dx' 

2ik J 


(3-81) 
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Take the average of Equation (3-81) 


<SP(x)> = 


N rL giklx-x'l 


N /•! 

•2: / 

J = 1 •'o 


<^f(x')5(x' -Xj))dx' (3-82) 


The mean of <'P(x') 6(x' - Xj)> is 


5(x' - Xj)> = J ^(x') 6(x' - Xj) f(Xj , x^, . . Xj^) dXj dx 2 • . . dXj^ (3-J 


where 

f(XpX2, . .Xj^3 = f(Xj) 

is the joint probability density function of the Xj. 

With the Xj statistically independent 

f(Xj , X2, . . . X^) = f(Xj ) fCXj) . . . f(Xj^ ) 


= n 

k= 1 

Equation (3-33) is therefore 


(3-84) 


<^(x') S(x' “ Xj)> = j ^(x') 5(x' - Xj) f(Xj^) dX|j (3-85) 

k= 1 


The product term can be expressed as 


N N 

f[ f(X{j)dxj^ = f(Xj)dXj Yl f(Xk)dXjj 
k = 1 k = 1 

k j 


Therefore, from Equation (3-85) 


/ N 

^(x') Y[ f(X]^)dxj^ 6(x' - Xj)f(Xj)dXj 
k = 1 


(3-86) 
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Note that 


N 

n f(xk)dxj, 

k = 1 

^ ^ j (3-88) 

= the field at x' averaged over all 
possible configurations of the 
otiier N-1 scatterers. 


<^(x') I x'> = II 


Then, from Equation (3-82) 
<^(x)> = e^^ - 


N gik|x“x'l r ^ 

or J J <^(x') I x'> 5(x' - Xj) f(Xj) dXj 


dx' 

(3-89) 


The integral over Xj is evaluated at the delta function, i, e. 

J mx')|x')5(x' - Xj)f(Xj)dXj = (❖(x')lx')f(x') 


(3-90) 


Thus, from Equation (3-89) 

N gikix-x'l 

<^(x)> = e^^ - o I <^f(x') jx')f(x')dx' (3-91) 

j = i o 

There are no j dependent terms in the integral, therefore the summation eval- 
uates the integral N times. 


('4'‘(x)> = 


- 



giklx -x'l 

<^(x') I x') f(x') dx' 

2ik 


(3-92) 


The rigorous solution of this equation requires solving the set of simultaneous 
linear algebraic equations for each scatterer and then performing the integrations 
to obtain <^(x') I x') , (see Eq. 3-88). This procedure is not feasible. 

The approximation introduced by Foldy consists of replacing <’5f(x') lx'), the 
field at the j th scatterer averaged over all possible configurations of the other 
N-1 scatterers, with <^(x')> , the average field at the jth scatterer position, 
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when that scatterer is not present. As Foldy notes, these averages differ only 
by a term of the order 1/N, thus if N is large, the approximation is reasonable. 

Therefore, with 


<^(x') I x'> ^ 


the mean field, from Equation (3-92), will be 





gik[x - x'l 

2ik 


f(x') dx' 


In addition, if the Xj are uniformly distributed, 

1 

f(x') = - 
L 


(3-93) 


(3-94) 


(3-95) 


Therefore 


(ijfpfx)) = 



giklx - x'l 

<’^(x')> dx' 

2ik 


(3-96) 


where A = N/L, 

The above integral equation represents the mean field for the slab under the 
Foldy approximation to the external field acting on each scatterer. 


(a) Series Expansion For the Mean Field and Coherent hitensity 


The Foldy Approximation to the mean field <^p(x)) , Equation (3-96), with the 
addition of the boundary conditions of the slab, represents the solution to a 
wave equation of the form 


_dx2 



ff X0?p(x)> 


or 


— - + k^d + aX) 
dx2 


<^Tp(x)> = 0 


(3-97) 


J ..i. .:l. 

A solution to Equation (3-9 1 ) is 

<^p(x)) = ^ ^ 

or 

<^p(x)> = (3-98) 

where 

K = k(l + (3-99) 

Thus, under the Poldy approximation for the mean field, the slab of discrete 
scatterers can be represented as a continuous medium with an 'effective' prop- 
agation constant K, given above. The mean field for x > L is then available 
from the transmission coefficient T, see Figure 3-3, 

(❖p(x)> = X > L (3-100) 

The mean field and coherent intensity (x)> (ijr* (x)) can nov/ be developed 
from the expansion solution for T to 0(c^ ), 

The transmission coefficient for a slab of width L and propagation constant K 
was derived in Appendix A, Equation (A-12), 

4k 



■ O L 

K = k{1 

Figure 3-3. Foldy Approximation Representation for the Mean Field 
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Expanding K to O(a^), (recall that cr may be complex), 


K = k(l + a\)'^ = k 


a\ (aX)^ 

1 + -h O(or^) 

2 8 


(3-102) 


Then 


ik^L -iki^ 

jiKL = gikL e 2 e ^ 


giKL _ gikL 


= ejJ'L 


oX (aXf 

1 + ik — L - k^ L2 

2 S 


1 - ik + ..0(c‘*) 


8 


(3-103) 


crX (ffX)2 (ffX)^ 

1 + ik — L - - ik 

2 8 8 


and 


gi2KL = gi2kL 


= pi2kL 


(nX)^ „ 

1 + ikaXL - 

2 

T (oX)^ - 

I + ikcrXL - k^ - ik 


4 ' 

1 - ik L...O(o^) 


4 


ioXf 


(3-104) 


Applying these expansions, we find 


(k + K)^ = k2(4 + 2aX - 


(o\f 


(k - K)2 - k^ 


iaXf 


L 4 J 


:>i2kL 


(k + K)2 - (k - K)2 ei2KL = k2 


j2 /„y ^2 

4 4 


ia\f (aX)^ , 

4 + 2aX 


4kKe*Kl' = 4k2 6^*^^ 


ffX cX (ffX)2 ^ (ffX)2 ^ (ffX)^ 

1 + ik — L + — + ik L - k^ 

2 2 8 8 8 . 
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The transmission coefficient is then 


•1 r Y ^ , X2 

jikL j — L-k^ — 

, V 2 2/ \ 8 ^ 

4k? 1 + ^ { "l6 " ~16 


This is of the form 


gikL I* j + Act + Ba^ ] 
[1 + aff + bcr^] 


where the A, B, a, and b are the coefficients in Equation (3-105), 


Expanding the denominator 


Therefore 


[1 + ao + =1 - an - btx^ + + O(ff^) 

= 1 - ao + (a^ - b)a^ 


T = + Ao + Bo^l [1 " aa + (a^ - b) 0^1 

= e^^[l + (A - a) o + (a^ - b - aA + B)a^] 


= 1 


\a (Xo)2 (Kaf 

1 + ik — L ik L 

2 16 8 


t2 t 2 i 

- k^ + e‘ 

8 16 


i2kL 


Rearranging terms 


... ko 

= 1 H- i — 

L 2 


k^ fl? 1 gizkL L 

XL X2(— + — ; 

4 \ 2 4k^ 4k^ 2ik 


The resulting mean field is therefore 




39 


or 


- e 



ka 

+ i — XL - 
2 




gi2kL 

4k^ 



(3-107) 


A comparison of the above result with the discrete solution expansion for the 
mean field, Equation (3-64), shows that the Foldy approximation represents the 
mean field exactly to O(a^), i,e. 


<^p(x)> = <^(x)> to 0(ff2) (3-108) 

The coherent intensity <’lrp(x}> <'^^=(x)> is found from Equation (3-107) and its 
complex conjugate, and the result will be identical to Equation (3-66), the dis- 
crete solution expansion. 

The Foldy approximation therefore represents the coherent intensity exactly to 
O(a^), i.e. 


|(^p(x)>p = K^(x)>P to 0(a2) (3-109) 

The total intensity expansion <^p(x) ^jt|(x)>, however, will not agree with the 
discrete solution expansion for (^f(x) 4f"*(x)) to O(a^). The Foldy approximation 
does not account for the incoherent intensity terra 

a CT* 

+ — XL 
4 

which is present in the discrete solution expansion, Equation (3-76), 

The implications of the Foldy approximation results developed here are dis- 
cussed in the following section. 


3-4. SUMMARY OF ONE -DIMENSIONAL MEDIA RESULTS 

The one-dimensional study of a slab of discrete scatterers presented in this 
chapter has produced a number of useful techniques and results which are ap- 
plicable to the volume scatter rain channel. They are summarized briefly below: 

(a) the Poisson distributed point scatterer, with an appropriate ’weighting 
function', can represent a distribution of discrete finite scatterers - 
the weighting function provides the physical characteristics of the 
scatterer for practical interpretation of results obtained. 
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(b) either discrete or continuous techniques may be employed for the de- 
velopment of the mean field and intensity expansions, as long as care 
is taken in the formulation to insure non- overlapping scatterers, 

(c) series expansion of the mean field to 0((T^) produces secular terms in 
L for both 0(a) and O(a^) - the slab width limitation for series expan- 
sion requires that 

2 

L < 

kloiX 

(d) the Foldy approximation expansion for the mean field agrees exactly, 
to O(a^), with the discrete Born series expansion. The Foldy approxi- 
mation can be used to represent the coherent intensity, however it is 
not applicable to the total intensity, where an incoherent intensity term 
of 0(a^ ) was found which is not accounted for in the Foldy expansion. 
The Foldy approximation has the same slab width limitation as given 
above. 

(e) an incoherent intensity term of 0{a ^ ) was found, which is secular in L, 
and which may be significant in the volume scatter case for certain 
wavelength/scatterer size ranges. 

The applicability of the one-dimensional results to the volume rain scatter case 
will be observed as the volume configuration solution is developed in the next 
chapter. 


CHAPTER 4 


PROPAGATION THROUGH A VOLUME DISTRIBUTION 
OF DISCRETE SCATTERERS 


A study of the multiple scattering of waves propagating through a volume of dis- 
crete scatterers is developed in this chapter. The scope of the study will be 
directed at determining a suitable representation applicable to the rain scatter 
channel and evaluating the statistical parameters of the field by discrete random 
media techniques, some of which were employed in the previous chapter for the 
one -dimensional case. 

The scattering volume is described by a uniform distribution of point scatterers 
whose position is described by a Poisson impulse process. The scattering char- 
acteristics of a single scatterer are represented, as in the one-dimensional ease, 
by an appropriate ^weighting factor'. Integral equations are developed for the 
mean field and intensity y the results are solved by Fourier transform techniques; 
and expansions for the far field are derived. The Foldy approximation to the 
mean field is applied and limitations on its range of validity determined. 


4-1. POINT SOURCE INCIDENT ON A VOLUME OF DISCRETE SCATTERERS 

We begin the development of the volume scattering study by defining the model 
used to describe the scattering medium. An isotropic spherical wave is assumed 
incident on a volume V of point scatterers whose position is determined by a 
Poisson impulse process. The isotropic incident wave was chosen because of its 
generality and ease of extension to other types of incident radiation such as the 
plane wave and beam wave. 

Consider the incident wave source point at , in an orthogonal coordinate system 
whose origin is located within the scattering volume V, as described in Figure 
4-1. The medium consists of N isotropic point scatterers randomly distributed 
in the volume V, with r^j the position vector of each scatterer (a random vari- 
able), defined by the Poisson impulse process 

N 

u(r) = ^ a5(r-rj^) (4-1) 

n = I 

where a is the complex weighting factor determined from the physical character- 
istics of a single scatterer, 

ttiJSGKWNG PAGS BL&m ^01 
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Figure 4-1. Coordinate Notation for Volmne Distribution of 

Discrete Scatterers 


The resultant field at the receiver point, ^(r), is desired. The field ^(r) obeys 
the reduced wave equation 


[V^ + n^(r)] 4f(r) = 0 


(4-2) 


where n(r) is the refractive index and k is the free space propagation constant, 
k == An time dependence is assumed. 


For the configuration defined above, the wave equation will be 

N 

E 

n - 1 




1 + a 5(r - r„ 


n) = 0 


(4-3) 


The above wave equation and Its associated boundary conditions convert to the 
integral equation 


^/(r,ro) ^ Go(r,rg) - k^ 


N r 
n = 1 •'v 


)6(r' - Tjj) W,ro)d?f' (4-4) 


where G(r, r') is the free space Green's function 

pikir - r'l 


G(r,r') = 


47Tlr-r'| 


(4-5) 
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and Gg (r, r^) is the incident spherical wave at i. e. 


Go(r, r„) 


gikir- rd 

47r!r - I 


The delta function evaluates the integral at r' = 


r^) = Gg(r,rg) - o 


E 


n = i 


G(r,r^)^(r^,ro) 


(4-6) 


(4-7) 


This equation represents the field at r as the sum of the incident wave in the 
absence of scatterers, and spherical waves, G(r, ), with an external field 
Tq) eminating from each of the N scatterers. 

The external field at the nth scatterer is the sum of the incident wave at Tjj plus 
the sum of the fields produced by all the other N-1 scatterers, i, e, 

❖(r„.r„) - Go(r„,ro) - G(rp,r^)GCr„,rp) (4-8) 

P = i 

p =/= n 


where G(tj^, tp) is the spherical wave produced at r^j from the scatterer at tp 
excited by the incident wave G(tp , ), 

Substituting ^(r^j , ) from Equation (4-8) into Equation (4-7) 


^(r,ro) = Go(r,ro) - a 



G(r, r„) 


N 




P - i 

p ¥= n 


or 


N 

4'(r,r^) = a G(r, r„) r^) 

n “ 1 

(4-9) 

N N 
n = 1 p = l 

where the prime in the p summation indicates "except for p = n," i, e. p n. 
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This result represents the fundamental multiple scattering equation for the field 
at r produced by a configuration of N Poisson distributed scatterers at in the 
volume V excited by a point source at . 

It can be shown that the resulting field, Equation (4-9), is also obtained by a 
Born expansion solution to 4f(r , ) where 


oa 

" S (4-10) 

j = 0 

and the expansion is carried out to O(of^)* This result was demonstrated for the 
one-dimensional case in Section 3-3(d); the procedure is similar for the three- 
dimensional case and need not be repeated here. 


Proceeding to the average of Equation (4-9) 

N 

E 

ill = 1 


mx, r^)) = Gjj(r, (7 / ^ G(r, G^(r^^, 


+ 


\n = 1 p = 1 / 


(4-11) 


The averaged summation terms above were developed in Appendix B. From 
Equation (B-25) 


2 **n) V = P f G(r, T^) Go(r„, r^) dr^, 

I n = 1 / •'v 


(4-12) 


and from Equation (B-33) 

N N 


iZ J^n) ^o) G(r„. tp^ 

n = I p = 1 


= ) f G(r, r^) G(tp, r^,) G(t^, v^) dr„ drp 

•'v '' (4-13) 


where 


N 

p = — 

V 


(4-14) 
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The mean field from Equation (4-11) is therefore 


<^r(r,r„)> = 


a p jf G(r, 


r') G(f , Fq) dr' 


k 4 J j G(r, r') G(r', r") G(r", r^) dr' dr" 


(4-15) 


(a) Mean Field Expansion by Application of Transverse Fourier Transforms 

The solution of the mean field integrals in Equation (4-15) are now considered 
for the case of a slab volume of width L in the direction of propagation z, and 
extending to ±00 in the x and y dimensions. The volume integrals can then be 
transformed into one-dimensional integrals in z by application of the two- 
dimensional Fourier transform in the transverse x - y plane. 

With 


r = + yay + za^ 


and 


let 


k = k^a^ + ky3y + k^a^ 


= xa^ + yay 


(4-1 6a) 


and 


The transverse Fourier transform pair for a function f(r) is then 


(4-1 6b) 


nz,k,) = jj 


. -ik* • r , 
f(r) e dr 


( 4 ^ 1 7 a) 


f(r) 




(4-1 7b) 
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where; 


dk^ = dk^dky,dr^ = dxdy 


and 


• r = kjj^x + kyY = k|. • 


For the Green's function 


G(r,r') 


_eik|r - r'| 
4irlr-!r'| 


the transfomi pair is of the form 


G(z 


-z’.k,)= ff 


G(r - r') e 


. -ikt * (r - r') 


dr* 


(4- 18a) 


G(r - r') = 


(IttY 


ff 


G(z - z', k^) ^ dkt (4-1 8b) 


The Green's function G(t, r') is a solution to the wave equation 

(V^ + k2)G(r,r') = 6(r - r') (4-19) 

ik * ii* 

If Equation (4-19) is multiplied by e ^ and integrated over dr^ to ±», the in- 
verse transform of G(r, r') is found to be (Tyras, 1969, p, 113) 

^ik 2 lz - z'l 

G(z - z', Iq) == (4-20) 

2ik, 

where 

- Ik^P (4-21) 

Utilizing the formulations of Equations (4-18) and (4-20), the Green’s functions 
in<^(r, Tq) , Equation (4-15), can be represented by their inverse Fourier 
transforms, i.e. 


<^r(r,r )> = 


1 


{27T)^ 

- or P 


// 


lk,lz-z„l 


2ik, 


Jhkfl 

s// 


ik, |z - z\ 


2ikz 


2il4 


dkj 


dr' 


+ 


"■'■// [s? // 

b// 


ik,|z - z'l 

® ik* • (r-r') 

e ' 

2ik, 


(4^22) 


dk* 


ik;|z'-z"| 


(27T)^ 


2ikl 


ik;'|z"-z. 


2ik! 




,ikr-(r"-rJdk;' 


dr dr'' 


The integral term over dr' in Equation (4-22) can be rearranged as 


Ij = -k^ G p 




^ik^lz - z'l 
2ikz 


dkt 


If 


Xlz'-Zcl -| 
2i!d ^ 


(4-23) 


J[kt • (i- - r ) + kt • (r - to)] 


>dr' 


With 


dr' = dr^ dz' 


the integral in Equation (4-23) over dr' is 
J"gi[kt • r - k; • rol ^-ir' • (kt - k^) ^ 


(4-24) 


Mt 


•r-k;.r„] J jz' 
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li , L.._ [ ..i.:: : 


Prom the delta function definition 


5(x-x')=— / 
(27T) 




we find that 


5(kt 


-k;)~ 

* (27T)2 JJ 


Therefore the bracketed integral term in Equation (4-24) is seen to be 


If 


dr; = (27 t) 2 S(kj - k;> 


Therefore Ij is 


Lk^lz-z'l ikl|z'-z^| 
e c 


- -k^ a p f 

(27T)2 J 


2ik^ 2ik; 


e‘'‘=*-r-'=;-'"J 6 (kj-k;)dk,dk;dz' 

The delta function evaluates one of the integrals at ~ k;, i.e, 

1 


I, - -k^ a p r 

^ (27T)2 


- ^ikjz-z'l gikjz'-z^l 

J 2 %. 2 ik^ 




or 


Ij ^ ~k^ a p f 

(27T)2 j 


ikz(|z - z i + |z -ZqO 


(27t)- 

^ik{ • (r-To) 


-4k| 


dkj dz' 


(4-25) 


(4-26) 


(4-27) 


(4-28) 


(4-29) 
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The integral term over dr' dr in Equation (4-22) can be rearranged as 




— ff 

C27T)^ J J 


Jk^\z - z'l 
2ik^ 


dkt 


// 


ik;|z'-z"' 


2iki 


dl4 


Xlz"-z,l 

2ik;' 


dkj, ^i[kt • r - k" . r J 


If 




(4~30) 


dz' 


The bracketed terms above over dri and drj' are seen, from Equation (4-26), to 
be delta functions in (k^ - k[) and (k^ - k^'). Therefore 


L = k^ p- 


5 ? // / 


ik^lz-z'l _ik;iz'-z"l ik;'lz"-z„| 

6 


7 I" "I 

e ^ e ^ 


2ik, 


2il4 


2K 


(4-31) 


^i[kt . r-k; • rj g^j^ _ 

The delta functions evaluate two of the dk^ integral pairs, at kj = and at 
k^' = k^ . First at k^' = k[ 


I2 = k^ p2 


If I 


ik,lz-z'l ik;iz'-z"l lk;|z"-zj 


(27t)2 j j j 2ik, 2il4 

gi[^ ^ *'ol 5^ - k|) dk^ dkj^ dz' dz 

then at 

ikj,(|z - z'l + |z' - z"l + Iz" - Zp I) 

r / f ' 

(27r) 

At * (r-^o) 


(4-32) 


I2 = k'^ or" p 


2 „2 


j?// 

dkj dz' dz' 


-8il4 


(4-33) 
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The results of Equations (4-29) and (4-33) can now be applied to Equation (4-22) 
to get the mean fileld 




■s? If 


! dk, 

2ik, ^ 


-k^ 


"5?// 

s// 


ik,(z - z 1+ |z -z 1) 


-4ki 


• (r - To) ^^ 34 ^ 


+ 


ikjdz -z 1+ \z -z 1+ |z -z^t) 

' ^ e“'*‘‘'-'"’dk,dz-dz’ 

(2tt)^ JJ -Sik^ ^ 


where 


= y/k^ - |k+l“ 


(4-35) 


O'" .Ider the second term integral over dk^ dz' in Equation (4-34) 




ik 2 (Iz - z'l + |z' - ZqI) 
e 

-4k| 


dz' 


(4-36) 


We select the origin, in rectangular coordinates, at the leading edge of the slab 
of width L, as in Figure 4-2. Then < z' < z, and the integral over z' above 
will be 



L 


^ikj,(z - Zq) 



dz' 




4 


k 


2 

Z 


Therefore 


-iff 


e^^z(^“Zo) ^ ^ 

^ ^ e'‘‘''f’^-’'“’dk, 


r r " ^o) "■ ~1 

J J 2ik2 ^ik^ J 


• (r - To) 


dk+ 


(4-37) 


(4-38) 
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Consider the third term integral over dk^ dz' dz'' in Equation (4-34) 


/ ik^Clz - 2 1 + Iz' - z 'I + |z' - ZqI) 

dz" (4-39) 




With z < z' , z" < z, 


(|z - z'l + Iz' - z"| 4 Iz" - z^l) = (Z - Zq) + (iz' - z"| - z' + z") 


and 



{r - To) 


ik2(z -Zq) 


-8il4 


rr 

O •' o 


dz" ^ ^ 

(4-40) 


The double integral over dz' dz" in Equation (4-40) is evaluated in Appendix C, 
Equation (C-2). The results are given by Equation (C-5). Then 


-/ 


jkz(z - Zo) p 


2ik, 


V- L 
+ 


u-84 


i8k| 


1614 


ik,2L"i 


16k^^ 




(4-41) 
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Finally, the resulting integrals andig, from Equations (4-38) and (4-41) are 
substituted back into Equation (4-34) to give 




“i//{ 


^ikjCz - Zg) 


2ik, 


1 - ap 


2ik, 


+ k" p2l 


L2 L 
+ 


ik,2Li 


^-Sk^ iSk^ 16kJ 16k 


(4-42) 


.,ikt • - r„) ^ 


From Equation (4- 16b) the term in brackets { . } above is seen to be the in- 

verse Fourier transform of r^)), i.e. (^(z-z^, kj)>. Let 


e 



p2 


4 


(4-43) 


Then 


^ik^(2 - zj 

(^(z - z kj)) = 

2ik^ 


k^ 

" L e 



e 


2 




— 

2ik| 


Le^ 




(4-44) 


Equations (4-42) and (4-44) represent the expansion solutions for the Fourier 
transform pair of the mean field of a semi-infinite slab of -width L, to O(e^). We 
now proceed to develop an asymptotic expression for (^(r, r^)) for large r, i.e. 
with the receiver in the far field of the incident and scattered waves. 


flD) Far Field Solution for the Mean Field 

Consider the integral expression for <’k(r, r^)) developed to O(c^) in Equation 
(4-42), We re-express the equation as 


<^(r, r^)) 


— ff 

(27t)2 J J 




(4-45) 


where 


B(kJ = 


ikz(z-Zo) 


2ik- 


1 - C7 p + 

2ikz 


L 


p2 „ — + 


1 


ik,2L 


) 


8 iSkj, 16k^ 16k| /J 

(4-46) 
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and 


= -s/k^ - Ik^P , dkj = dkj^ dky 


The field quantities axe now represented by an^ar displacement parameters, 
as defined by Figure 4-3. From Figure 4-3a 


and from Figure 4-3b 


Note also that 


> +0O 


(x - Xq) = R cos Q 

(4-47a) 

(y ^ yo) = R sin 0 

(4-47b) 

kj. = k^ cos 0 

(4-48a) 

ky = kj. sin 0 

(4-4Sb) 

dfcj^ dky = k|. dk^ d0 

(4-48c) 

* +CO /* °* C 

(4-49) 


dy 


kjdkj 




-7T 


Therefore 

kt • (r - r^) = k^(x - x^) + ky(y - y^) 

= cos 0 R cos 6 + kj sin 0 R sin 6 
= Rk^[cos0cos0 + sin 0 sin 0] 

= Rk^ cos (0 - 6) 

The integral over dkj in Equation (4-45), called I, is then 

-+ir 


I 


"/ I 

J O •f -SI 


d0 B(k^) e 


i R cos (0-0) 


(4-50) 


(4-51) 




With , B(k 2 ) is re-ejipresBed as 


- k^2 fz _ z^) 


BOc^) = 


2iVk^ -kf 


k^a2p2 / i2 


1 + i 


ff p L 


2^/k^^ 

^ , i2L \/k2 -k? 

L 1 ~e ^ 


'(k^-kj)^ 8 ^ sVk^ -kf 16(k2-kJ) yj 

Equation (4-51) is rearranged as 


/•«> r , 

1 = 1 k, B(k,) j e’^ ‘ 
J 0 -n 


cos (0-0) 


d0 


dk+ 


The term in brackets is of the form of the zero-th order Bessel function 
(Abramomtz, 1964), 


J„{z) 


1 . 

= -/ 


z COS |3 


d<3 


or 


1 . 


z cos j3 


dj3 


Let 


cos (0 - 0) = cos |3 , |3 = 0 - 0 , d|3 

The integral limits are 

0 — ■Hr”^j3 = 7r-0 
Since the integrand is an even, periodic function 


= d0 


/ +7T ^ 7T - 0 /* +JT 

d0 ^ / dj3 I dp 

•7T -IT - fl -IT 


(4-52) 


(4-53) 


(4-54) 


(4-55) 
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Therefore the term in brackets of I, Equation (4-53), is 





i R k* cos (^ - 6) j . 
e ^ d(p 



e* P dp 


= (27r)JoCktR) 


(4-56) 


and 


. = 2ir f ktB(kt)Jo(kt R)dkt 
o 


(4-57) 


Substituting I above back into Equation (4-45) 


<4r(r,r )> = — /* 
27t7„ 


kt B(k^) Jo(kt R) dkt 


(4-58) 


This integral can be converted to one extending from to +<3° by expressing the 
Bessel function in terms of Hanliel functions (Tyras, 1969) 


Therefore 


Jj^(z) = i [H^^^(z) + H^^>(z)] 


(4-59) 

(4-60) 


J^(ktR) =-[H^‘^(kt R) + R)1 


(4-61) 


and 


= -H[,‘>(kt R) 


(4-62) 


Since k^Ris real, the H[,^^(ktR e-i^^) term is equivalent to the term 

along the negative real axis, 0 to -«». Therefore, from Equation (4-58) 


i /•+“ 

<'k(r, r^))=— kt B(kt) >(1^ R) dk^ 

47T ^ _oo 


(4-63) 
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The integral of Equation (4-63) is evaluated asymptotically in Appendix D for 
large r by; (a) expressing the Hankel function in its asymptotic approximation, 
(b) employing appropriate coordinate transformations, (c) evaluation of saddle 
points, and (d) determination of a steepest descent (constant phase) path of 
integration. 


Application of the results of Appendix D, Equation (D-36), to Equation (4-63) 
gives the far field expansion for the mean field, 


aiki r 


<’^(r)> ^ 


4itrr 


1 + i 


kff p L 
2 cos a 


8 cos^ a 


ko^ L 
8 cos^ a 


Idcos'^ck! 16 cos'^ a 


(4-64) 


where a is the an^e between the axis and the receiver point r, as seen in 
Figure 4-3. 

The far field expansion contains one first order term, and secular terms in L 
and L A comparison of the far field expansion above with the one-dimensional 
result, Equation (3-64), shows that the terms are similar, and with a = 0, the 
forward scatter direction, the terms are identical in form. Note also that for 
L = 0, 


gikr 

<^(r)> = (4-65) 

4iin' 


which is the far field incident wave approximatloa. 

(c) Total Intensity Expansion 
The total intensity 

<|^(r, r„) |2> = <4r.{r. r„) ^*(r, r^)) (4-66) 

is now evaluated, to O(n^), for the integral field representation developed in 
Section 4-1. From Equation (4-9) 


N 

4r(r, = G^(x, Tq) - a ^ G(r,rjj) G^^(r^, r^) 

n = 1 

N N 
n - 1 p = 1 


59 


i.. 


the oomplex conjugate is found 

= G«(r.r„) - a* ^ G*(r, r„) GJ(r„. r„) 


N 


n = 1 


(4-67) 


N N , 

Y, E’ G*(r.r„)G*(rp,r„)G»(r„.rp) 
n = 1 p = 1 

The product is 

N 

'Kr,r„)4'»(r,r„) = G„GJ - k^oG* G(r. r„) G„(r„, r„) 


n = i 


N N 


H-k-lcr^G* £ £'G(r,r„)G(rp,r„)G(r„,rp) 

n = 1 p = I 


N 

I 

n = 1 


-k2a*G„ £ G»(r,r„)G5(r„,r„) 


(4-68) 


N N 


+ k^(< 72 )«G„ £ £ G*(r,r„)G»(rp,r„)G*(r„.rp) 

n = 1 p = 1 

N N 

+ k^ac7* 12 G*(r, r„0 G*(rjj, r^) 


n = 1 


n - 1 


Taking the expectation value of Equation (4-68) 

I 


(nr,r^)^*ir,r^)) = <t'k(r, r„)P> = 


(47t)2 ir-r^P 
N 


- ff G^' / Y, G(r. r^) G^(v^, r^) 


(4-69) 


,D = 1 


+ G* 


/ N N ^ \ 

\12 12 ^fr’'*n3G(rp,ro)G(r„.rp)\ 

\n = 1 p = I / 
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N 




,n = 1 


N N / 

z: E 

n = 1 p = 1 
N N 


+ k‘'(a2)*G„( ■£ J]’ G*(t.r„)G*frp.r„)Gnr„.rp)> 


+ k*aa-*(2 G(r.r„)G„0-„,r„) G*(t. r^) GJ(rp. r„)^ 

The eajjectatioti values of the summation terms were evaluated in Appendix B, 
The first and third (single) summations are found from Equation (B-25), the 
second and fourth (double) summations from Equation (B-33), and the last (un- 
primed) summation from Equation (B-31). The resulting total intensity is then 


<l^(r,ro)P> = 


(47T)2lr-rQ!^ 


oG^p f 

J \f 


- k^aG^p I G(r,r')G„(r'.ro)dr' 


- k^ or* Gq p / G*(r, r') G*(r', r^) dr' 

\r 


+ k^ G* p2 y J G(r, r') G(r", r^) G(r', r") dr' dr" 

+ J j G*(r, r') G*(r", r^) G*(r', r") dr' dr" 

a*[p^y yC(r,r' 


(4-70) 




) G^(r', r^) G*(r, r") G*(r", r^) dr' dr" 


+ p r G(r,r')G„(r',rJG*(r,r')G*(r',r„)dr' 
*^v 
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The solution of Equation (4-70) is accomplished by application of the transverse 
Fourier transform to each term, as was done for <4r(r, r^)) in Section 4-1 (a). 
With each Green's function replaced by its transform from Equation (4-18b) 

C(r. r') = yj G(z - z', k,) e‘‘‘> ‘ dk, 

and evaluated, a solution of the form of Equation (4-42) is obtained. A compari- 
son of Equation (4-70) with Equation (4-15) shows that the first two single integral 
terms for the intensity are identical in form to the single integral term for the 
mean, and the first two double integral terms are likewise the same as the double 
integral mean term. Thus the results of the mean field integral evaluations can 
be directly applied to solve the first five terms of Equation (4-70). Since the far 
field solution is desired for this study, the explicit development of the intensity 
expansion for the far field will be considered. 

Before solving for the far field total intensity, however, the coherent intensity 
will be developed and compared with the results given above, and the last inte- 
gral in Equation (4-70) will be seen to represent an incoherent intensity term. 

(d) Coherent Intensity 

The coherent intensity 


K^rCr.r^))!^ = (^(r.r^,))<^nr,ro)> 
is found from Equation (4-15) 

<^(r, r^)) = Gq(f, Tq) - a p r G(r, r') G(r', r^) dr' 

+ CF^ y* J G(r, r') G(r', r") G(r", r^) dr' dr" 
Taking the complex conjugate 

(4F*''(r, r^)) = G*(r, r^) - k^ 0 ='’ p f G*(r, r') G*(r, r^) dr' 

Jv 

+ k'^(a^)* J' J G*(r, r') G*(r', r") G*(r", r^) dr' dr" 


(4-71) 


(4-72) 
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The coherent intensity is the product 


imr, ro))|2 = <4r{r, r^)) r^)) = G* - a p G* f G(r, r') G(r', r^) dr' 

*^v 

-k^n^'pG^ f Gnr,r')G=^(r'.r^)dr' 

J V 


+ k'*' CT^ G* U G(r, r') G(r', r") G(r", r^) dr' dr" (4-73) 

+ k'‘(n2)=i= p2 Gq J* G*(r, r') G=<'(r', r") G'^Cr", r„) dr' dr" 


+ 


ao^p^ J J G(r, r') G(r', r„) G'^Cr, r") G=*'(r", r„) dr' dr' 


Comparison of this result with Equation (4-70) for the total intensity shows that 
all terms are identical except for the' last term in the total intensity, which is 
not present in the coherent expansion. Therefore, to O(a^) 

(mr, ro)l^> = U^(r, r^,)>l2 + k^ a p / G(r, r') G^(r', r^,) G'"(r, r') G*(r', r^,) dr' 

•'v (4-74) 

The last term of Equation (4-74) is of interest, since it represents an incoherent 
intensity term, ,, in the total intensity expansion. will be evaluated in 
the next section by asymptotic expansion for large r. 

(e) Incoherent Intensity 

The incoherent intensity, from Equation (4-74), is 


line = k" a 


a*p f G(r.r')G(r’,r„)G»(r,r")G*(r",ro)dr' (4-75) 

J \r 


with r' - r". 


63 



We proceed to solve by application of the transverse Fourier transform to 
Equation (4-75), as was accomplished in Sectica 4-1 (a) for the mean field. Re- 
placing each Green's function in Equation (4-75) by its inverse Fourier trans- 
form, Equation (4-lSb) 


C 7 ff p /* r rPT, / , X ikt • (r - r') „ 

^nc - ~ j8 jj tG(z - z,kt)e t dk^ 

G(z" - z„,k;")e‘‘‘‘ ' dr' 

The inverse functions are replaced by their one-dimensional equivalent solu- 
tions, from Equation (4-20) 


(4-76) 


^nc 


k^ CT a* p f ~ ^-i^z '!z” - Zq! 


(27T) 


8 


/ 


2ik, 


2ikl 


-2ik; 


-2ik: 


"2 ““’■Z -“'■2 -■ “"Z 

i[k^ • (r - r') + k; . (r' - r^) + k”(r - r") + k;" • (r" - rj] 


(4-77) 


dr' dk^ dk; dk[' dk|' 


The k^ exponential can be expressed as 


-ir' . (kt - k[) -ir" • (Ic[' - k'/') m, + k;') . r - (k; + kl") . r J 


(4-78) 


Since r' = r" 



[(ktTk;')-(k;+kni 


dr^ = (27T)2 5i(kj + k;') - (k; + k;")] 


(4-79) 


Therefore 


^inc 


k4 


ff CT* p r 

27t1^ J 


i[k2[z~z l + kjiz - ZqI - k^ |2 - z [ - |z -z^l] 


(27T) 

^i[(kt + k;') • r-(k; + k;") • Tq] 


I6k,k; k"i4" 


(4-80) 

6(k^ - k;) 5(k|' - k;") dk^ dk[ dk^' dk^" dz' 
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The delta functions evaluate the k^. integrals at and k^'= k[" 


^inc " 


k4 


a ff'" p I e 


/ 


i[k2(z - z') + k^^lz' -Zq! ■“ kjlz -z'l -k^'lz' - 2^|] 


(27T)4 
i[(kt + kt') • (r-To)] 


2 i,»2 


16 k. 


dkj dfc;' dz' 


With Zg < z' < z 


(4-81) 


^nc 


C27T)'’ 


/ 


16k3 k"^ 


z * z 


k^ aa*p C 

-Zo) 

2 i,n2 


^l[(k,tk").(r-r,)l 

(4-82) 


,4 .. f. ik_(z-z„) -ik"(z-Zo) 

k^cF a^’p / e ° e ^ ° ;r/T, +k"'»*fr-r 

\.C - j "^^4 / 1^' 


(2irr J IS k| kj 

Evaluating the integral over dz' 


Jo 


( 4 - 83 ) 


a a'^‘ p f Zo)^ik 2 (z z^,) +k'‘ *(r-r )1 

I.„ = L j (r rj] 


(27T) 

Rearranging terms 


16k^-k;'2 


(4-84) 


^nc 


k'^aa=>=p 


1 6(27t)'^ 




^ik,.(r-t„)d^ 


-ik^(z - Zo) 


k;- 


^ik;-(r-Odk; 

(4-85) 


The two integrals in Equation (4-85) are seen to be transverse Fourier trans- 
forms of conjugate functions, i. e, (see Eq. 4- 17b) 


/ 


jk^(z - Zq) 


^ikt • Cr-ro)(jj^ = f(y _ 


(4-86a) 


and 


/ 


-ik^(z - Zq) ^ 

2 gik; • (r - r„) ^ ^2^)2 fi:(r „ r^) 


(4-S6b) 
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Therefore, from Equation (4-85) 


L 


me 


a a* p L 
16 


f(r - r^,) F(r - r^) 


(4-87) 


Since 


Then 


f(r - r^)f*(r - r^) = |f(r - r^)!- 


^inc “ 


k** a a* p L 
16 


If(r - ro)P 


(4-88) 


The function f (r - ) is now evaluated in the far field by the method of stationary 

phase. With 


re - to) “ 


ikS! 


- ^o) 


(4-89) 


The above integral is seen to be identical to the integral in Equation (4-45), with 

ik2(z - z„) 


B(kJ = 


’7 ■' T 


(4-90) 


which was evaluated in the far field in Section 4-1 (b), and Appendix D. 

The method is applied here to f(r - r^,), with only the major steps listed. From 
Equation (D-4) 


4 -+0O 


ik,R 


Vk;B(u.)c-dk. 


(4-91) 


where and R are as defined in Figure 4-3 . With 


k^ = k sin W 


k^ = k cos W 


\/i^ = Vk^rTw" 
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Therefore 


ik*R 
e ' 


B(kt) = 

dk^ = 


gikR sin W = gikr sin W sin a 
gikr sin W sin ft 
cos^ W 
k cos W dW 


,-‘f 


f(r - r„) = 


%/S 7F^ R 


/ 


Vic %/ sin W 

k^ cos^ W 


k cos W dW 


1 _iiL /'V^iirWe^*^>^‘=‘'®Ov -a) 


(4-92) 


Stt^ k R 


/ 


cos W 


dW 


The SDP determined in Appendix D is applicable for the above integral also, and 
the asymptotic representation given by Equation (D~27) can be applied. With 


qCWg) = k, W3 = o!, lq"(z)| = k 


f(W,) = 


VsIrTa 
cos a 


Therefore, from Equation (D-27), for large r 
f(r - r^) ^ 


1 /ZffVsinoi irk + ii 

e / — e 


Stt^ k R 

dkr 


rk cos a 


2 7T k r cc a 


and 

lf(r - = [4 7 t’ k^ cos^ 

Then, from Equation (4-88) 


(4-93) 


(4-94) 


me 


k^ a a* p L 
64 k^ cos^ a 
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■with 


era* = (a^ + a?) 

_ k2(aj +ff?)pL 
64 7T^ -p‘ cos^ a 


(4-95) 


(f ) Far Field Solution for Total Intensity 


The total intensity (l^(r, was found, to 0(a^}, in Section (d) as, Equation 

(4-74) 


where is the incoherent intensity term given by Equation (4-75) and by 
Equation (4-95) for large r. 


The far field approximation for the total intensity can now be developed from the 
far field expansion to the mean field, Equation (4-64). The complex conjugate 
of Equation (4-64) is 


,-ikr 


<^*(r)> 


-4i7T r 


kff*pL k^(a^)*p^L^ k(a“)*p~L 

I _ I ;; + j 


(o^)* p^ 


2 cos O' 

n2 


8 cos^ oc 


8 cos^ a 


(4-96) 


+ ^ ^ cos a 


1 6 cos*^ or 16 cos'^ a 


The coherent intensity, to O(a^) is then 


<4f(r)) <^*(r)> 


I 

1 6 7T^ 


k a p L 

1 + i 

2 cos a 


k^ 

8 cos^ a 


k a^ p^ L 

8 cos^ a 1 6 cos'^ a 




7 


1 6 cos*^ a 


gi2kL cos tt 


k a* p L 
2 cos « 


k^ a (j* p^ k^(a^)*p^ k(a^)*p^ L 

+ - + i 


4 cos^ a 8 cos^ a 

(g2)#p2 (a-)*p^ 


8 cos^ a 


,-i2kL cos ot 


16 cos^ « 16 cos'^ar 
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Or 




1 


16 Tp- 


1 + i 


+ — ( a a’^' - 

4 cos^ a \ 

(a^ + (cj^)*) 

16 cos'^ a 


k{ff - a*) p L 
2 cos a 

+ (a^ )*’\ k(a^ - (ct^ )*) p^ L 


2 

n2 


- 1 


8 CQS^ 


a 


(4-97) 


1 6 cos'^ (X 


gi2kL cos cc ^ cos a) 


Rearranging terms 


<^(r)) = 


1 


1 6 7T^ 


k C; p L k^ cr? p^ 1? 

1 ! + 1_ 

cos a 2 005“^ Di 


k o. Ut p^ L p^ ff. O: 

+ — I — L + ^ sin (2 k L cos a) 


2 cos^ a 1 6 cos'^ a 

(oj -(7?)P^ ^ P^C(T? -fff) 
8 cos'^ a 8 cos'*’ a 


(4-98) 


cos (2 k L cos a) 


The addition of the far field result for from Equation (4-95), results in the 
total intensity far field expansion 


(l^(r)|2) = 


1 


16 


k o, p L k^(o^ +cr?) p L 
1 1 ^ ^ 


cos a 


4 cos^ o: 


k^ of \J' k ff_ Oi p^ L m Oj p^ 

+ + — + — sin (2 k L cos a) 

2 cos^ a 2 cos^ a 16 cos'^ o: (4-99) 


(aj-opp^ {af-af)p'^ 


cos (2 k L cos a) 


8 cos*^ « 8 cos*^ a 

For Of = 0, i. e. the f^'Tward scattering direction. Equation (4-99) gives 

<l^(r)|2> r - * , 




{ \r 

- -1 1 - k p Oj L + k^ p^ a? — + 


k^ p(ff j + a?) L k p^ Oj L 

— + — — 

4 J 2 

(4-100) 


0.0: p^ (of - of)p^ 

+ -1-5 — sin (2 k L) (1 - cos 2 k L) 


16 


8 
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where the term in brackets [ 3 ^ incoherent intensiiy term. 


The far field total intensity expansion contains one first order term, and secular 
terms in L and L^. A comparison of Equation (4-100) and Equation (3-76), the 
one-dimensional result, shows that all terms of the expansion are identical in 
form. 

Also, for L = 0, 


<l^(r)^> = 

16 

which is seen from Equation (4-65) to be the far field intensity of the incident 
wave. 

The applicability of the Eoldy approximation to the volume scattering case will 
now be investigated. 


4-2. FOLDY APPROXIMATION FOR A VOLUME OF DISCRETE 
SCATTERERS 

The Foldy approximation to the external field acting on a volume of discrete 
scatterers is applied here to the volume slab configuration described in the pre- 
vious section. The mean field r^)) is developed, and the resulting ex- 

pansions compared to the mean and intensity expansions derived in the previous 
section. The development is similar to the one- dimensional slab case presented 
in Section 3-3. 

Consider the wave equation for the spherical wave incident on a volume of dis- 
crete scatterers, as described in Figure 4-1, and by Equation (4-3) 

|v2 + k2 

The wave equation and its associated boundary conditions convert to the integral 
equation 


N 

\ + a Y, SC*’ - *‘n> 


n = I 


} 


^(r, r^) = 0 


(4-101) 


N f 

>k(r,ro) = Gg(r,r^) - k2 a J G(r, r') 5(r'-rj,) ^(4 r^) dr' (4-102) 

n = 1 V 


where 


Go(r-ro) = 


-e 


ik|r- r^l 


47r|r-rQ| 


and 


G(r.r') 


•\ - 


-giklr - r'l 
47T|r-r'| 


The average of Equation (4-102) Is 


<^(r,r^)> = G^(r,r^) - a ^ f GCr,r')(5(r' - r^) ^(r', r„)> dr' (4-103) 

n = l*' V 

The mean of S(r' - r^) N[f(r' , r^) is 

<5(r' - Fj^) «f(r', r^,)) = j 5(r' - t^) 4f(r', Tq) f(r^ , .... ) dr^ , dr^ .... dr 

J fArr 


(4-1 04) 


where 


f(ri,r 2 ,...rj^) = f(r„) 
is the joint probability function of the . 

With the Fjj statistically independent 


N 


= f(r,)f(rj)...f(r„) = f[ f(fp) 

P= 1 


(4-105) 
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Therefore Equation (4-104) is 


<6(r' - r„)^(r\r„)) = f 6(r' - [] (4-106) 

p= 1 

The product term can be expressed as 


N 


p= 1 p = 1 

p n 


(4-107) 


Therefore Equation (4-106) is 


N 

4r(r', r^) [] f(Fp) dFp 5(r' - ) f(r^) dr^ 

P = 1 

p =?^ n (4-108) 


<W,r^)5(r' - r„)> = jjj 


Note that 


// / 4r<r'.r„) f(ip)drp = 


P 

P n 


(4-109) 


which is the field at r' averaged over all possible configurations of the other 
N-1 scatterers. Therefore, from Equation (4-108) 


<^(r',r^)5(r' - r^)) 



<W,rQ)lr'>6(r' - 


The integral is evaluated by the delta function at r' = , i.e. 

{nr\r^)5{T' - r^)) = (4^(r', ro)|r') f(r') 
Therefore from Equation (4-103) 


(4-110) 


(4-111) 


<4^(r, r^)) = Gg(r,rQ) - a 



G(r,r'X^(r',ro)ir^>f(r') dr' 


(4-1 12) 
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■t : 


d:.--... -I. .. ■ .J 


r 


There are no n dependent terms in the integral, so the integral is evaluated N 
times by the summation 

<^(r,r^)> = G^(r,r^) - k^uN 

The approximation introduced by Foldy consists of replacing <^(r', the 

field at the nth scatterer averaged over all possible configurations of the other 
N-1 scatterers, with r^)) , the average field at the nth scatterer position, 

when that scatterer is not present. These averages differ only by a term of the 
order l/N, thus for large N, the approximation is reasonable. 

Therefore with 




G(r, 


t^)\r‘)f(r') dr' 


( 4-113 


the mean field, from Equation (4-113) will be 


( 4 - 114 ) 


<^Tp (r, r^)) = G 3 (r, r^) - u N T G(r, r') <4^(r', r^,)> f(r') dr' 

J v 


Also, since the scatterers are uniformly distributed in V, 


1 

f(r') = “ 
V 


( 4 - 115 ) 


Therefore 


where 


<^p(r,ro)> = GQ(r,ro) - k^ a p 



G(r,r') mr',ro))dr' 


( 4 - 116 ) 


P 


N 

V 


The above integral equation represents the mean field for the volume V of dis- 
crete scatterers under the Foldy approximation to the external field acting on 
each scatterer. 
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(a) Mean Meld Expansion 


The Foldy approximation integral equation for the mean field (r, r^)) , Equa- 
tion (4-116), mth the addition of the boundary conditions of the scattering volume, 
is seen to represent a wave equation for <^p(r, r^)) of the form 

[V2 + k“ ] <^p(r, r^)> = -k2 a p r^)) ' (4-117) 


or 


+ k2(i + ap)] = 0 (4-118) 


The above equation defines an 'effective^ propagation constant, K', for which the 
mean wave propagates in the scattering medium. The mean wave, with an in- 
cident spherical wave, is thus represented in the slab as 


<^p(r.ro)> 


^iK'ir-rol 


4TT|r-rQl 


(4-1 19) 


where 


Expanding K' to O(a^) 


K' = k(l + a pf^ 


or 


K' = k 


ap (ap)2 

+ — 

2 8 „ 


K' 


ap (ap)^ 

k + k k 

2 8 


(4-120) 


(4-121) 


The results thus far are seen to parallel directly with the one-dimensional case, 
with K as given in Equation (3-102), and X replaced by p. 


Under the Foldy approximation for the mean field, the volume slab of discrete 
scatterers can be replaced by a 'continuous* medium with an effective propagation 
constant K', as displayed in Mgure 4-4. The field expansion for ('^pir, r^)) in 
Region III, z > L, is desired. 
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!■ ■■ ■■ ■■ V - 4 I _ I .-.a....... ; --i. 


X 



K' = k(1 +op)'^^ 

Figure 4-4. Foldy ApproximatLon Representation for Mean Field 
of a Volume Slab of Discrete Scatterers. 


The solution for <^p(r, r^^)) is reduced to a one -dimensional evaluation by con- 
sidering the transverse Fourier transform in each region, i. e. 


• 


~ e ‘ “ _iV 

kj)> = + r e ' 

2ik, 

z < 0 

(4- 122a) 

* 



. 0 < z < L 

(4-1 22b) 




, z > L 

(4- 122c) 


where 

= y/k^ - 


(4-123a) 


k; = x/k' 2 - kj (4-123b) 
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The boundary conditions at z = 0 and z - L are 


<^J> = <^„> , z = 0 

(4-1 24a) 

<^jj> -= , z = L 

(4-124b) 

9 ^ 9 

^ - 0 

(4- 124c) 

9 r«J 9 ^ 

— ^ 2 = L 

9z “ 9z 

(4-124d) 

The resulting equations from the boundary conditions are 


A + B - r - 

(4-125a) 

= 0 

(4-12 5b) 

g-i^z^o 

AK;-BK; + rk = 

-Z z z 2i 

(4-1 25c) 

AK^giK^zL _ - Tk^e^''^ = 0 

(4-125d) 

For ease of notation, as before. Equation (4-43), let 


ap (np)2 

e = — = 

2 4 

(4-126) 

Then 


K' = k + ke-k — 

(4-127) 


2 


and 


k; = _ 1^2 = + 2 e k2 - k2 


= \/k| + 2ek^ = k^ 



(4-128) 
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Expanding the squaj?e root term to O(e^) 




1 + 


e-l(^ 

Ak 


= gis 


where g is the term in brackets of Equation (4- 129a). 
The expansion for the exponentials in is 


iKiL 


ik_ 

= (a ^ 


t • ft)’ ■ ■ i (t) 


(4-129a) 
(4-1 29b) 


and 


= e 


ikjL 


Jr2 1 

1 + i — e L r {k~ + i L) 

k, 2 k! 


= pe 


ik^L 


(4-130a) 


(4- 130b) 


-iKiL 
e * 



L 


C4-131) 


where p is the term in brackets, Equation (4-13 Oa), and p* denotes its com- 
plex conjugate. 

The expansions of K^, and are applied to the four simultaneous 

equations, Equations (4-125) to give 


^-ik2Zo 


A + B - r 

(4- 132a) 

2ik2 


pA + B - T = 0 

(4-13 2b) 

g-ikz^o 


gA - gB 4 r = 

(4-132C) 

2ik^ 


gpA - gp^’' B - T - 0 

(4-132d) 


The solution for the transmission coefficient T was developed by a Cramer’s 
reduction of the determinent for the four simultaneous equations. Equations 
(4-132). 


The result, in series form to O(e^), was found, and the resulting field 
from Equation (4'-122c) is available. The inverse transform of the field, in 
Region III, from Equation (4-122c) and the T ejqpansion is 




,ikz(z “ Zo) 


2ik, 


1 , 

1 + i — Le-i Le^ 

L k, 2 k| 


(4-133) 


1 


1/k 


1/k 


r + — [ — 1 e 

2 k| 4\kJ 


i2k,L 0 


A comparison of the above result with the discrete solution expansion for the 
mean field inverse Fourier transform, Equation (4-44), shows the two expans- 
ions to be exactly equivalent. Thus the Foldy ^proximation expansion to the 
mean field for the volume scatter case has been shown, as in the one-dimensional 
case, to represent the mean field exactly, i. e. 


<4fp(r,rQ)> = <^(r,r^j)) to 


(4-134) 


The far field expansion for the mean developed in Section 4-1 ^), Equation (4-64), 
is therefore also applicable to the Foldy expansion result. 


(b) Intensity 

The results of the previous section have shown that the mean field <S^(r, r^)) can 
be represented, to O(n^) by <^p(r, r^^)), thus the coherent Intensity (^^Cr, r^)) 
r^)) will likewise be represented by the FoMy approximation, to 0 ( 0 ^). 

While the Foldy approximation represents the coherent intensity, it does not, 
however, account for the Incoherent intensity term, found in the far field ex- 
pansion, Equation (4-95), as 

(Oj + o? ) p L 

64 7T^ cos^ a. 


which is present in the total intensity far field expansion. Equation (4-99). 

The Foldy approximation to the mean field is therefore found to be a useful tech- 
nique in the representation of the mean field and coherent intensity’ by consider- 
ing the volume slab of discrete scatterers as a continuous random medium with 
an ’effective’ propagation constant. However, the Foldy approximation does not 
represent the total intensity, where an incoherent intensity term of O(n^) was 
found which was not accounted for in the Foldy expansion. 
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CHAPTER 5 



.a-. . . . 4 


I 


IMPLICATIONS OF DISCRETE SCATTERING RESULTS 
TO MULTIPLE SCATTERING IN RAIN 


The two previous chapters have investigated the intensity of a wave propagating 
through a distribution of discrete scatterers by application of random media 
techniques. In this chapter, the study concludes with a discussion of some im- 
plications of the results obtained to the rain scatter channel. 

First, some general comments can be made on the mean field and intensity ex- 
pansions. The analyses of Chapter 3 have demonstrated that either discrete or 
continuous techniques may be employed for the development of the mean field 
and intensity expansions, as long as care is tal<en in the formulations to Insure 
non-overlapping scatterers. The resulting expansions for the mean field were 
shown to be equivalent to a first order Born series expansion solution, and the 
limit of convergence for the expansions is equivalent to the Born expansion 
limit. 

The mean field and intensity expansions for a distribution of discrete scatterers, 
developed in Chapters 3 and 4, have shown that the volume distribution results 
proceed directly from the one-dimensional distribution findings. The similarity 
of the volume scattering results to the one-dimensional case was not unexpected, 
however the term by term agreement, even to the incoherent intensity term, in- 
dicates that a one-dimensional evaluation could provide useful quantative as well 
as qualitative predictions of multiple scattering effects in rain scattering. 

Next, the applicability of the isotropic scatterer assumption to the rain scatter 
problem Is discussed. The volume distribution results of Chapter 4 were ob- 
tained for isotropic point scatterers, and extension to the rain scatter case is 
dependent on the validity of this assumption for the case of a distribution of di- 
electric spheres. The induced or scattered field produced by a dielectric sphere 
in the presence of an incident electric field can be shown by a direct solution of 
spherical wave functions of the incident and induced waves to be a vector field 
with components which depend on the angle of observation to the sphere (Stratton, 
p. 566), Hence the spherical rain drop cannot, in general, be considered as an 
'isotropic scatterer'. 

If, however, the discussion of the rain scatterer is limbed to the Rayleigh scat- 
tering region and to the forward direction only, the scattered field can be shown 
to b6 a scalar function of the incident field (Stratton, p. 435), and the results of 
Chapter 4 can be considered. 
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Thus* the question of how the multiple scattering results developed in this study 
relate tc the classical 'single scattering' method can be considered by comparing 
expansions of both methods in the forward scattering direction, where angular 
dependence is not a factor. 

A summary of the classical technique used to determine rain attenuation is pre- 
sented in Appendix E. The second order series expansion for the classical 
'single scattering' intensity of a wave propagating through rain is given in Equa- 
tion (E-14}, When this expansion is compared with the forward direction ( a = 

0) multiple scattering expansion given by Equation (4-100), the 'single scattering' 
formulation is seen to be similar in form to the first three terms of the multiple 
scattering expansion. All other terms, including the incoherent intensity term, 
are not included. 

It is also of interest to note that the Foldy approximation expansion was found 
to represent all terms of the coherent intensity (Section 4- 2b), thus the Foldy 
approximation appears to be more representative of the total intensity than is 
the classical 'single scattering' formulation. 

It must be emphasized again that these observations are made for the forward 
scattering direction only, and the complete vector scatter solution may produce 
additional factors which are not evident in the scalar results discussed here. 

In conclusion, the study has demonstrated the feasibility of using discrete random 
media techniques for the determination of multiple scattering effects in propa- 
gation through a volume of discrete scatterers, and has provided some insight 
into the more general problem of multiple scattering in a volume of rain. 
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APPENDIX A 


TRANSMISSION COEFFICIENT FOR A ONE- DIMENSIONAL 
RANDOM MEDIUM SLAB 


Consider a slab of width, L, with a propagation constant K, with an incident wave 
■*1^0 (^) as shown in Figure A-1. In regions I and III the propagation constant is k, 
and transmission coefficients F, A, B, and T, as shown in the figure, are 
assumed. 

The field in region I is 

^j(x) = + Y^^-ikx (A-1) 

In region II 

^d(x) = Ae^^ + (A-2) 


k K k 



0 L 

Region Region Region 


I II III 

Figure A-1. One-Dimensional Slab Configuration 


A-1 




In region III 




The boundary conditions at x = O and x = L require that 

^i(P) = ^n(O) 

^ii(L) = (L) 


(L) 


From the boundary condition equations 


Rearranging terms 


1 + r = A + B 


Ae>KL gg-iKL = 

ik - ikr = iKA - iKB 

iKAe’KL _ iKBe-'KL = 


A + B - r = I 
Ae*^'- + Be~'*^L = T 
KA - KB + kr = k 

KAe'KL _ K3e-iKL = {jy 


Multiply Equation (A-4) by k, and add to Equation (A-6) 

(k -h K) A + (k - K) B = 2k 


A-2 


Therefore 


2k k - K 
k+K“k+K® 


(A-8) 


Multiply Equation (A-5) by k, and subtract Equation (A-7) 

(k - K) AeiKL + (k + K) Be-iKL = o 


Therefore 


A 


k + K 
k - I< 


g-i2KL g 


(A-9) 


Equate Equations (A-8) and (A-9) 


2k(k-K) 

(k-K)2 _ (k + K)2 e-i2KL 


(A-10) 


Substitute B into Equation (A-9) to find A 

-2(K + k)ke-'2*^L 

A 

(k-K)2 - (k+K)2 e*2KL 


(A-il) 


Finally, substitute Equations (A-10) and (A-11) into Equation (A-5) to find T, the 
desired transmission coefficient 

^ ^ -4kKe‘^^- 

(k-K)2 - (k + K)2 e"i2KL 


To express T as a positive function, multiply by -e'^KL / _gi2KL ^ 


4kKe^^^ ^ 

(k + K)2 - (k-K)^ e‘2KL 


(A-12) 
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APPENDIX B 


EVALUATION OF MEAN FUNCTIONS FOE ONE-DIMENSIONAL 
AND VOLUME DISTRIBUTIONS OF DISCRETE SCATTERERS 


In the development of the mean fields and intensities produced by a distribution 
of point scatterers, the evaluation of means for a series of random functions are 
required. This Appendix develops the means for these functions, for the one- 
dimensional case. Section (a), and the volume distribution, Section (b)» 

(a) One-Dimensional Distribution of Scatterers 

Consider h(Xj), a function of the random variable set Xj. The mean value of 
h(Xj) is; 


<h(xp) = J' h(xj) f(Xj , X2 , . . . Xj^ ) dxp dx2 , 


dx 


N 


(B-1) 


where f(xi, X 2 , . . . x^) is the joint probability density function of the . If 
the xj are statistically independent and uniformly distiibuted along O to L, 


f(xi, X 2 , . . . x^) == fCxi)f(x2) fCx^) 

= 11 1= fiY 

L L ” ' ' L VL/ 


(B-2) 


Consider the function 


N 

n ° I 


(B-3) 


From above 




f(xi, Xj, X3, . . . x^) dxi dx2 . . dx^ 


(B-4) 



Ih(xi) + h(x2) + . . . + h(x^i )1 dxj dx2 . 


. .dx„ 


Each term of the series is of the form 


|h(x,,dx, //-/.n 


dxs 


(B-5) 


j 4= i 


There are N of these terms, therefore 




N 


n dXj (B-6) 
j = 1 

i 


Since 


n L ^ 

... f n<‘>‘j = LN-‘ 

•'O j = 1 


(B-7) 


j ^ i 


Therefore 




h(xp dxj 


N 

L 



h(Xj) dxj 


(B-8) 


Let X = N/L, the uniform number density then 


2 h(x„)\ = X I h(Xj) dXj 

n = 1 / Jo 


(B-9) 


Consider now a function h(xi , X2) of two random variables of the set xj , 

<h(xi,X 2 )) = I h(xi ,Xj ) f(Xi ,X2 , . . . X^ ) dxi dx2 . . . dx^j 


■iWI 


(B-10) 


h(Xi,X2)dxj dx2 . . . dxN 


B-2 


Similarly 







<h(Xi,X2) + h(x3,X4)> = ( “) / [h(xi,X2) + h(x3,X4)l dxj dx2 . . dx^ (B-11) 


and 



“ N N 

,E E liCxj, xj) 

j = 1 1=1 


dx^-^ dx2 . • . dXj^ 


Jh(x^,Xj)dXj. jj j n 


m =r= j 


+ (N^ - N) J‘ h(Xj , Xj) dxj dxj // /, n 


m j, I 


i. e, , the first term is for N terms where j = 1, second term is for terms 

where j f 1. 

Since 


and 


n L ^ L N 

J •' 0 m = 1 


m j 


(B-13) 




f[ 

m = 1 
m ^ j, 1 


(B-14) 


Therefore 



E 


1= i 



. 


N(L^ 


■>/ 


h(3<j, Xj) 


dX; 


+ (N2 





h(Xj , Xj) dXj dxj 


(B-15) 


N 

L 



+ 


-N 



h(Xj, Xj) dXj dXj 


For large N 


N2 - N 
L2 



= \2 


(B-16) 


Therefore 



N \ 

1 = 1 / 



(B-17) 


The first integral results from the N^-N terms of the series where j 4 = 1, while 
the second integral results from the N terms where j = 1, i. e. , where the 
scatterers are allowed to overlap. 


B-4 


For the case where scatterers are not allowed to overlap, only the first integral 
is present, 1, e. , 


\ j = 1 1 = 1 / -^0^0 

where the prime ' denotes j r 1. , 

(b) Volume Distribution of Scatterers 

In this section, the same summed averages as above are derived for the volume 
distribution case. The development is nearly identical, and only major results 
are presented. 

Consider, h{ Rj^ ), a function of the random variable set . With the R^^ 
statistically indep endent and uniformly distributed in the volume V, 

f(ri,r2,...rN) = fdi ) Rr.) . . • f(rN ) (B-19) 

and 

i=l,2,...N (B-19a) 

Then the mean value of h(Ri) is evaluated as 


<h(R] )}= I h(rj ) f(ti ) Kr, ) . . . fCr^ ) drj drj . . . dr^ 

Jy 

(B-20) 

J h(ri ) drj drj . . . drj^ 



and 


(h(Ri) + h(R2)> = 



[h(r|) + h(r2>l dtj dr2 . . . dr^, 


N 


th(rj) + h(r2) + . . . + h(rN)] dr^ dtj. . .dr^j 




[Ii(ri)dri] dt2 dtg . , . 


+ 


J [h(T2)dr2] drj dr3 


dri 


N 


Each term of the series is of the form 

Jh(ii)dxi jj- J n 

j ^ i 

There are N of these terms, therefore 

(t«v) -( 7 )” »/■«•.//■ -/.n-i 

i 


Since 


Therefore 




j = 1 
j^i 


drj = V^-i 



(B-21) 


(B-22) 


(B-23) 


Let the number density be p i, e. , 


P = 


N 

V 


(B-24) 


then 


n = 1 / -^V 


(B-25) 


Consider now a function h(R^, Rj ) of two random variables of the set , 


(hCRj.Rj)) 


= /hfri. 





.V 


1^2) P(J *1 ) PCr^) . * . p(%) drj dr^ . _ drj^ 


h(rj , r^) drj dr^ . * . dr^^ 


(B-26) 


Similarly 







(hlRpR^) + h(R3,R4» = (-1 j [h(Ri,R2) + h(R3,R4)l dr^ dr2 . . . dr^^ 

(B-27) 

The sum is therefore 


N N \ 

E E «Rn'V) 

= i p = i / 


■©7 


N N 


E E ’’(%.rp) 

Ln = 1 p = 1 


dti dr^ . . . drj^ 




+ N J h(r^,rp)drj, iHh 


j n J 


j n,p 

(B-28) 


B-7 


The first integral consists of the (N^ - N) terms where n 5 ^ p, the second 
integral the N terms where n = p. 


Since 

ff -Jn 

j n,p 

and 

//-■■/ n. 

j ^ n 

Therefore 


(B-29) 


(B-30) 






h(tn, tp) dr„ dtp 


+ p 



h(r„,r„) dr„ 


(B-31) 


where, for large N 


-N ^ Ni = p 2 

y2 y2 


(B-32) 


For the case where the scatterers are not allowed to overlap, i.e. j 1 , only 
the first integral is present, i.e., 



N / 


-I — 1 


li(R„,Rp) 



h(rn, Fp) dTj, 


dr^ 


(B-33) 


where the ' indicates n 7 ^ p. 
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APPE3SIDIX C 


SOLUTION OF DOUBLE INTEGKALS REQUIRED 
FOR FIELD SOLUTIONS 


Two integrals which, occur in the analysis of the mean field and coherent intensity 
of a slab of width L will be evaluated here. 

Consider the first integral, 


«/ 0 o 


gik(Ix - X I 4- |x " X I + X ^ x) 
for X > L, X > x' , x", see Figure C-la 


r r 

J Q J O 




gik(|x'^x"|-i-x"-x') dx'dx" 


CC-I) 


(C-2) 


The integral is evaluated in two regions, I and II, as shown in Figure C-lb. k 
Region I, x' > x" and (lx' - x"l + x" - x') = 0 and the integral limits are 



(C-3) 


In Region II, > x' and ( lx' - x" I + x" - x') = 2x” - 2x' and the integration 
limits are 



(C-4) 


Therefore: 


/•L ' 

/"L 

f dx' 

x' + j dx" 

F 

o 

A' 


C~1 




Continuing the evaluation 


fL 

pik2L _ pik2x'“ 

I dx' 

x' + e-i^ 2 x' 

/o 

2ik J 


= x' dx' + f 

J 0 Jo 


h gik2L e-ik2L _ j 

2ik 


dx' 



j2 ^ik2L 

.1. 

>ik2L _ f 


L 

2 

2ik 1 

_ -2ik _ 


2ik 

V- 

1 

J. 

gik2L 1 

4 


2 

4R2 

4k^ 2ik 



The second integral to be evaluated 


(C-5) 


I 


2 



g-ik(lx - x'l + lx' - x"l + x” - x) dx" 


CC-6) 


is seen to be the complex conjugate of Ij . The evaluation therefore follows 
directly from Equation (C-5). 


I? J_ 

2 4k^ 


,-ik2L 


4k^ 


L 

+ — 
2ik 


(C-7) 
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ASYMPTOTIC EVALUATION OF MEAN FIELD INTEGRAL 


Consider the integral given by Equation (4-63), 

1 

h = — I (D-1) 


and the required radiation condition 

ImVk^ -kj > 0 (D-2) 

To evaluate for large r, we apply the asymptotic s^proximation for large argu- 
ments for (Abramowitz #9. 2. 3) 


H^'CkjR) 



(D-3) 


with 


-V < aig(lq R) < 27T 


therefore, Equation p-1) is 


where 


= \/k^ - k| 
or 

k^ = yk^ - kf (D-5) 

The square root function kj; is a multi-valued function of k^, and its properties 
are investigated in a two plane ''Reimann sheet*' diagram for kj., shown in Fig- 
ure D-l(a). The top sheet, characterized by Im k^ > 0, is the proper sheet, 
where the radiation condition is maintained; the bottom sheet, where Im kg; < 0, 


D-1 



is an improper sheet. In order to cross from one sheet of the Rieman surface 
to another, the cut, called a branch cut, must be crossed. 

In the complex planes, the sheets -will have branch cuts along Im = 0, and 
branch points at k^ = ±k (k real). 

The path of integration P is shown on the top sheet of Figure D-l(a) where an 
infinitesimal amoimt of loss has been introduced (Im k > 0 ), which results in a 
slight displacement of the branch points at ±k. The four quadrants of the two 
kt sheets have been labeled Tj , . . . T 4 , . B 4 . 

To eliminate the branch points, we introduce a new complex variable, W = 

Wj + i Wj, and define a W-plane by the transformation 

= k sin W (D-6) 

= k(sin Wj, cosh Wj -t- i cos Wj. sinh Wp (D-S) 

and also 

k| = ~ k^ = cos^ W (D-9) 

k2=kcosW (D-10) 

where the + root is chosen to allow k^ = k when = 0. The kj. plane is thus 
mapped into the W-plane as shown in Figure D-1^). k^ = k sin W is single 
valued, and the path of integration P proceeds as shown, with each sector of the 
W-plane representing the various quadrants of the k^-plane as labled. 

We replace the displacement (z ~ Zq) with its angular equivalent, (see Fig. 4-3), 

(z - Zq) = rcosa (D-11) 


and 


R = rsin a CD-32) 

Therefore, Equation (D-4) in the W-plane can be developed 

\/k^ - - k-s/l - sin^ W 


“ k cos W 


■s/ic^ - y/k sin W 


^iktR _ gin _ ^ijtr gjji w sin a 


dkj. = k cos W dW 
And, from Equation (4-52) 

I 


BCI4) 


- pikr cos W cos a 


a pL 


2ik cos W 4 cos^ W 


2 


- 1 


ka^ p 


iP' 


I _ gi2Lk cos W 


(D-13) 


- 1 


2 cos^ W \ 8 8 k cos W 16 k^ cos^ W 


or 


B(kp = TCcos W) “ (D-I4) 

where T(cos W) is the bracketed term [ ] in Equation (D-13). Therefore, 

Equation (D-4) is 


L = 


e 

Vs 


3 4 ^+« 

■' / Vk^Vsin W T(cos W) ^ * k cos W dW 

(D-15 : 


Let 


Then 


F(W) = VsinW T(cos W) cos W 


(D-16) 


I 


1 



F(W) dW 


F(W) dW 


(D-17) 


The above integral can be evaluated for lart'e r asymptotically by deforming the 
integration path P into a path of steepest descent, SDP, through the saddle 
points. 


D-4 


To determine the SDP, we need to examine the phase characteristLcs of the ex- 
ponent in Equation (D-17), since the steepest descent path is a constant phase 
path in the neighborhood of the saddle points. Recall that, (see Eig. 4-3) 


iri = V(x - + (y - + (z - z^Y 


(D-18) 


r = lrle^“, < a < — 


(D-19) 


The exponent of Equation (D-19) is then 


ikr cos(W - a) = ik — + e ” “3 

2 


— + i w j - a) ^ g-i(Wr + i Wi - a)] 


(D-20) 


2 


! 

Re[ikrcos(W - a)] = k — - e sin(W- - ce) (D-21) 

2 


Im[ikrcos(W - a)] = k — - e "^0 cos(W_ - a) 

2 ^ 


(D-22) 


As seen on the W-plane plot, the integration path begins at Wj = +“ and ends at 
Wj = -°®. We require that the integrand of , (Eq. (D-17)) decays at tliese end 
points, i, e, , 


Re[ikrcos(W - a)] ->• 0 as Wj ^ ±°° 


From Equation (D-21), this is found to require 


7T 7i 

w -5- < ft < — as Wf ->■ +o® 

^2 2 ^ 


(D-23a) 


V—: -t: 


(D-23b) 




7T 


7T 


W_-» <a< — asW. 

"2 2 ^ 


Since -?r/2 < a < tt/ 2, the above conditions are seen to be met and the integrand 
will decay at the end points of the path defined on the W-plane plot, since the 
path remains within the allowable region. 

The steepest descent path, SDP, is a constant phase path, i.e. , 

Im[ikrcos(W - a)] = constant as Wj “5- ±oo 
From Equation (D-22), this requires that 

7T 

W_ -5- a , as W; ^ 

2 ' 

ir 

W- a H — , as Wi -oo (D-25b) 

1 « 1 


(D-24) 


(D-25a) 


We see that the ends of the path lie within the allowable regions of the W-plane 
since ■ jt/ 2 < a < tt/2 by its definition (Fig. 4-3). Thus the path P can be de- 
formed into SDP, and the integral can be evaluated asymptotically. Figuie D-2 
shows the Pgj^ where a small amount of loss a « 4 - i e has been introduced to 
move the poles at ±Wp to allow P to be defoi’med into P^^ without introducing 
a singularity. 

The method outlined by Felson and Marcuvitz 1973 (Chapter 4) will be used to 
develop the first order approximation for Ij , (Eq. (D-17)). Felson and Marcuvitz 
(Eq. (4. 2a), (lb)) show that an SDP integral of the form 


im = 


can be asyniptotically represented by 


X 

PSD 


f(z) dz 


(D-26) 


m ) - 



f(z„)e 


iSl^(zs) 


if 


as oo where Zg is a first order saddle point, (real) and 

q(z) = i^(z), ^real 


(D-27) 
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For Equation (D-17), Q = r 


q(W) = ik cos(W - a) 
^(W) = k cos(W - ct) 


(D-28) 


(D-29) 


The saddle point, W^, is found at q' (W) = 0 

q'(W) = -ik sin(W - a) 


W, = « 


CD-30) 


Therefore 

q(Wg) = k 

q '(W) = -k sin (W - a) 
^'(Ws) = 0 

q"(W) = -kcosCW - a) 
q"(Ws) = -k 


(D-31) 

(D-32) 


f 

(D-33) 


f(W) = F(W) = v^ihTWTCcosW)cosW 
f(Wg) = Vsin a 


1 (7 p L k 

— + + j — 

2ik 4 cos a 1 6 cos^ a 


(D-34) 


L 

4 — + i 


ff2p2 




- 1 


j2kL COS ce 


1 6 cos^ a 32 k cos'^ ce 32 k cos^ ce 


The asymptotic representation for Ij, Equation (D-27) is then found from Equa- 
tion (D-17) with the above substitutions; 



— Vsin a T(cos ce) 
rk 


jkr + i j 


- (D-35) 


where the + sign for the iff/4 term is selected to provide zero argument at r = 0, 
and T(cos a) is the bracketed term [ ] in Equation (D-34). Recalling that 


D-8 




D-9 



APPENDIX E 

SUMMARY OF CLASSICAL METHODS IN RAIN 
ATTENUATION AT MILLIMETER WAVELENGTHS 



APPENDIX E 








SUMMARY OP CLASSICAL METHODS IN RAIN 
ATTENUATION AT MILLIMETER WAVELENGTHS 


We begin the development of the rain attenuation discussion by briefly reviewing 
the classical techniques used to predict rain effects on electromagnetic propa^ 
gation. Most previous studies of rain attenuation, including the often referenced 
works of Gunn and East (1954), Medhurst (1965), Setzer (1970), Hogg (1968, 1969), 
Oguchi (1960, 1964), and Hogg and Chu (1975), among many others, have used the 
so-called 'single scattering' approximation either directly or through the assump- 
tions made in their analyses. 

While the method does not rigorously postulate single scattering or specifically 
exclude multiple scattering effects, the assumed exponential decay for the wave 
intensity, the use of cross-sections based on forward scattering only, and the 
assumption of Independent drop statistics with superposition of intensity, all 
lead to a model which does not account directly for multiple scattering. 

The classical 'single scattering' formulation assumes an exponential decay for 
the intensity of the wave as it propagates through a rain volume. Let 4^ ( r ) be 
the wave propagated through the rain, and 4'o (r) the wave in the absence of rain. 
The intensity of the wave is then of the form 


1% (r)l^ 


(E-1) 


where p is the drop density, the number of drops per unit volume; is the 
attenuation cross section of a single spherical drop, defined below; and L is 
the path length. 

The attenuation cross section of the single spherical raindrop is defined as 


Qt = 


Total power extracted from the incident wave (watts) 
Incident power density (watts/m^ ) 


where the units of Qt • 


E-1 


I 


i 


- 1 


I. 


The cross section of a spherical drop can be expressed in terms of the Mie co- 
efficients and (Stratton, p. 569) 


27T \ 

^77 Zj ^ I)Re[an + bn] 


k2 


n = 1 


(E-2) 


where the an and bn are complex functions with arguments dependent on drop 
radius, wavelength, and m, the complex refractive index of the water drop. Qt 
is temperature dependent through m, which varies by as much as a factor of 
100% over the temperature range of 0 to 40'’C (Ray, 1972), 

For a distribution of drops with varying radii, 


pQt = 


yQt(a)n(a)da 


(E-3) 


where n(a)da is the number of drops per unit volume with radii between a and 
a + da, and the integration is over the range of drop radii. 

The drop size distribution, n(a) , has been determined by direct observation to 
be approximately proportional to the rainfall rate, R, by an exponential relation 
of the form (Laws and Parsons, 1943; Marshall and Palmer, 1948) 

n(a) = n(o)e-«aR-'2I (E-4) 

where n(0) and « are constants determined from fits of measured distributions. 

The resulting intensity for the case of a distribution of drop sizes, at the wave- 
length X, for a rainfall rate R is then, from Equation (E-1) 

l^(r) P _ -L f Qt(a, X) n(a, r) da 

■ — e ^il:— ol 

l%(r)|2 


If the rainfall rate varies along the path L, the path variable must also be 
integrated, i.e. 


i’FfrlP _ - dx y*da Qt(a, X) n[a, R(x)] 

!% (r)P 


(E-6) 
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The intensity ratio above is usually expressed as a dB ratio, referred to as the 
’attenuation' of the wave. From Equation (E-6) 


A(dB) = iOlog 


10 


^da Qt (a, X) n[a, R(x)] 


(E-7) 


= 4.343 y* dx yQ^(a,X)n[a, R(x)] 


da 


The classical or 'single scattering' formulation for rain induced attenuation was 
shown to be of the form 


Iv^(r) P ^ ^-pQ^L 


where the total cross section, in terms of the Mie coefficients, is 


CO 

Qt =-T z] + DRela^ + b^l 

n = 1 


Consider the region where 2 ttsl « X , i, e. , the Rayleigh scattering region, 
where the drop size is much smaller than the incident wavelength. The Mie 
coefficient expansions then converge very rapidly, and only the first two co- 
efficients need be considered. As developed by Stratton (1941, p. 571), 

a. (E-8) 


2 . - 1 /27raV ^ 1 . - 1 ^27ra V 

3 + 2 \ X y 15 + 2 \ X y 


(E-9) 


When the ratio of radius to wavelength is so small that (2 7ra /X)^ can be neglected 
with respect to (2 tt a /X)^ , only bi is significent. Then the cross section will be 


lir 2ir 

Qt = — 3Re(bi) = — 3 
k2 k2 


-(ka)3 Im ( — ~ 

3 -i- 2/ 


(E-10) 
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where 



Therefore 


Qt 


= k4ira3 Im 


“ 1 
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^ + 2 , 


= kft 


where 


ft 


47ra3 Im 
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2 _ 1 


,m 
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The intensity ratio is then 


I ^(r) _ g-kp^jL 

I (r) P 


With Ik p)3i L I <1, which corresponds to attenuation values of 4.343dB or less, 
the eicponential can be expressed in series form, i.e. , 


I q'Cr) P 

I % (r) 1^ 


kpjSjL + 


k2p2^2L2 

2 


+ 0([p^jL]3) 


(E-14) 
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